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An  approximate  set  of  invariant  functions  for  the  dressed  vertex  amplitude  was 
found.  An  asymptotic  solution  to  the  unrenormalized  Schwinger- Dyson  equations  of 
Quantum  Electrodynamics  was  obtained  which  joined  smoothly  with  the  solutions 
found  by  a  perturbation  technique.  The  photon  propagator  is  approximated  by  its 
form  near  the  mass  shell.  The  vertex  equation  was  separated  from  higher  order 
members  of  the  hierarchy  at  the  second  order  in  the  coupling  constant  with  the  aid 
of  H.  S.  Green's  generalization  of  Ward's  Identity.  No  infinities  were  substracted 
to  obtain  the  solutions.  The  function  multiplying  the  matrix  7'^  is  found  to  be 
dominant  everywhere. 


CHAPTER  I 
INTRODUCTION 

Before  the  arrival  of  quantum  mechanics,  there  were  two  major  problems  that 
seem  to  point  out  the  existence  of  flaws  in  the  classical  theory  of  electromagnetic 
fields.  They  were  the  blackbody  radiation  and  the  theory  of  electrons  in  atoms.  The 
introduction  of  the  finite  quantum  of  action  by  Planck  successfully  resolved  the  first 
of  these  problems.  The  difficulties  connected  with  describing  the  electron  and  other 
entities  as  "particles"  are  much  more  serious.  Basically,  the  nature  of  the  treatment 
of  the  electron  as  a  classical  particle  with  finite  mass  and  extent  leads  to  a  dilemma. 
If  we  assume  the  electron  to  be  a  point  without  structure,  then,  in  classical  theory, 
the  total  energy  of  its  electromagnetic  field  becomes  infinite  implying  an  infinitely 
massive  particle.  On  the  other  hand,  if  we  assume  the  electron  to  have  a  finite 
extent,  its  interaction  with  the  electromagnetic  field  created  by  its  own  charge 
distribution  will  produce  stresses  tending  to  explode  the  charge. 

The  quantum  theory  of  the  electron  and  the  electromagnetic  field  was  intro- 
duced in  the  1920s  by  Heisenberg  and  Pauli,  Dirac  and  others.  Unfortunately,  they 
were  not  able  to  solve  the  equation  that  they  derived  in  a  manner  that  was  free  of 
infinite  quantities.  Although  the  energy  of  the  electromagnetic  field  diverged  as  the 
logarithm  of  the  electron  radius  in  the  quantum  theory  instead  of  as  its  reciprocal 
as  it  does  in  the  classical  theory,  a  second  infinity  appeared  that  is  associated  with 
the  charge.  This  infinity  appeared  in  integrals  that  diverged  quadratically  as  the 
upper  limit  increased. -"^ 


According  to  Dirac's  hole  theory,^  the  creation  of  an  electron-positron  pair 
by  a  photon  may  be  interpreted  in  the  following  way.  The  vacuum  consists  of 
an  infinite  sea  of  negative-energy  electrons.  A  photon  may  be  absorbed  by  this 
vacuum  raising  a  negative  energy  electron  to  a  positive  energy  state  and  therefore 
creating  a  hole  in  the  vacuum.  The  latter  will  appear  as  a  positron,  so  that  we 
obtain  an  electron-positron  pair.  Thus  the  vacuum  may  be  considered  as  a  sort  of 
"polarizable"  medium,  because  it  potentially  contains  electron-positron  pairs.  A 
photon  may  now  interact  with  this  polarizable  vacuum  even  if  its  energy  is  not 
sufficient  to  create  a  real  electron-positron  pair.  In  this  case  only  a  "virtual"  pair 
is  created  and  this  annihilates  soon  afterwards.  Although  the  apperance  of  virtual 
pairs  improved  the  divergence  of  the  self-energy,  it  introduced  new  problems  that  do 
not  have  their  counterpart  in  the  classical  theory.  With  the  introduction  of  virtual 
pairs,  it  was  found  that  the  polarizability  of  the  vacuum  is  infinite. 

The  Dirac  theory  also  predicts  that,  for  hydrogen-like  atoms,  states  with  the 
same  total  quantum  number  n  and  angular  momentum  j  are  degenerate  (same 
energy).  It  was  noted,  however,  that  the  polarization  of  the  vacuum  discussed  in  the 
preceding  paragraph  would  split  this  degeneracy.  In  particular  the  2^5i/2  and  the 
2^Pi/2  levels  should  be  separated  by  a  small  amount.  In  1947,  Lamb  and  Retherford^ 
made  a  direct  measurement  of  this  separation.  They  found  that  the  2'^5'i/2  is  above 
the  2'^Pi/2  by  1058  megahertz.  Actual  theoretical  calculations  using  a  form  of  a 
perturbation  expansion  of  the  equations  of  quantum  electrodynamics  (QED)  of 
the  splitting  gave  rise  to  divergent  integrals.  Bethe^  circumvented  the  problem  by 
simply  limiting  the  range  of  integration  over  the  divergent  integrals.  Bethe  reasoned 
that  at  energies  larger  than  the  rest  energy  of  the  electron,  relativistic  effects  become 
very  important  and  have  to  be  included.  These  inclusions,  at  the  end,  will  amount 
to  the  introduction  of  a  cut-off  in  the  integration  limits.  The  final  results  were, 
therefore,  dependent  on  a  cut-off  parameter.   This  dependency  was  logarithmic  in 


the  cut-ofF  parameter  and  thus  insensitive  to  the  actual  value  of  the  parameter. 
This  parameter  could  be  described  as  the  maximum  energy  of  a  photon  that  was 
emitted  and  absorbed  by  the  atomic  electron  of  hydrogen.  He  set  this  maximum 
energy  to  be  the  rest  energy  of  the  electron.  With  this  technique,  Bethe  arrive  at 
approximately  the  value  measured  by  Lamb  and  Retherford  (the  so-called  "Lamb 
shift"). 

The  fundamental  equations  of  QED  can  be  written  in  different  forms.  One 
way  of  doing  this  is  by  writing  the  perturbation  expansion  in  a  series  of  powers  of 
the  electron  charge.  This  expansion  can  be  written  in  terms  of  propagators  (Green's 
function)  as  developed  by  Feynman^  and  interaction  sites  at  which  three  "particles" 
can  interact.  Schwinger^  and  Dyson^  derived  an  infinite  hierarchy  of  integral  equa- 
tions that  describes  the  interaction  of  an  electron  with  the  radiation  field.  In  the 
future,  we  will  refer  to  this  hierarchy  as  the  Schwinger-Dyson  hierarchy  as  it  is 
usually  known.  The  equations  in  this  hierarchy  can  be  solved  by  straightforward 
iteration  to  obtain  the  perturbation  series.  This  perturbative  approach  always  leads 
to  infinite  integrals. 

Numerous  attempts  were  made  to  eliminate  the  divergencies  in  a  rigorous 
manner  in  the  period  following  the  invention  of  quantum  field  theory  around  1925 
until  after  World  War  II.  This  problem  was  solved  about  the  time  of  the  Lamb 
and  Retherford  experiments,  by  the  development  of  the  renormahzation  theory. 
The  idea  of  renormahzation  can  be  interpreted  as  a  rearrangement  of  the  pertur- 
bation series  so  that  the  new  series  converges  and  has  finite  terms.  This  rearrange- 
ment amounts  to  expanding  the  electron  and  photon  propagators  and  their  join- 
ing vertex  around  the  mass  shell  (i.e.  for  values  of  the  photon  energy-momentum 
near  zero  and  electron's  square  energy-momentum  near  its  experimental  square 
mass).  The  infinity  associated  with  the  charge  is  removed  by  rescaling  the  prop- 
agators, wave  functions  and  vertex  parts^'^  with  the  introduction  of  the  so-called 


renormalization  constants  Zi,Z2  and  Z3.  They  are  associated  with  the  vertex 
function,  electron  propagator  and  photon  propagator  respectively.  Attempts  to 
calculate  these  constants  using  perturbation  theory  have  led  to  the  conclusion  that 
these  constants  must  be  infinite!  The  renormalization  theory  emphasized  the  covari- 
ant  aspects  more  strongly.  Schwinger,  Tomonaga^°  and  Feynman,^^  independently 
developed  the  first  Lorentz  covariant  schemes  designed  to  eliminate  the  divergencies 
in  a  more  acceptable  manner  (the  renormalization  theory). 

Renormalization  theory  has  enjoyed  a  remarkable  success  in  the  calculation 
of  numerous  effects  such  as  the  Lamb  shift,  the  anomalous  magnetic  moment,  the 
hyperfine  structure  of  the  hydrogen  atom,  and  other  relativistic  phenomena.  Quan- 
tum electrodynamics  has  become  a  model  for  other  field  theories  to  follow.  It  is 
important  then,  to  study  the  underlying  mathematical  structure  of  QED  in  or- 
der to  better  understand  field  theories  in  general.  The  success  of  the  non-Abelian 
gauge  theories  in  unifying  the  electromagnetic  interaction  with  the  weak  interaction 
further  encourages  the  efforts  to  understand  and  resolve  the  ambiguities  of  QED. 
The  current  theories  of  electro- weak  and  strong  interactions  are  based  on  the  same 
underlying  mathematical  structure. 

If  the  mathematical  techniques  (perturbation  theories)  used  in  QED  were  com- 
plete and  satisfactory  theory,  it  would  be  a  mathematically  rigorous  and  logically 
consistent  structure  which  allows  at  least  in  principle  the  calculation  of  all  radiative 
processes.  Because  calculations  have  been  based  on  the  rearrangement  of  infinite 
series  that  may  not  remove  all  of  the  infinite  quantities  (.Zi,  Z2,  and  Z3  for  exam- 
ple), this  cannot  be  said  without  reservations  of  the  theory  in  the  present  state  of 
development. 

There  are  two  major  points  of  view  regarding  the  infinities  found  in  QED. 
One  states  that  there  is  something  fundamentally  wrong  in  the  foundations  of  the 


theory.  The  other  assumes  that  all  the  difficulties  arise  from  the  use  of  inadequate 
mathematical  methods  in  solving  the  fundamental  equations  of  the  theory.  The  first 
point  of  view  states  that  the  solution  to  the  problem  must  be  found  in  a  modified 
theory  or  in  a  completely  new  one  (if  possible).  The  arguments  against  this  radical 
approach  are  found  in  the  actual  success  of  renormalization  theory  in  calculating 
such  effects  as  the  Lamb  shifts,  the  anomalous  magnetic  moment  of  the  electron, 
and  the  fine-structure  constant.  In  order  to  get  a  flavor  of  the  kind  of  agreement 
that  can  be  accomplished  with  this  theory  the  following  results  are  presented:^'^ 

1.  Lamb  shift  in  hydrogen 

8E^^p  =  1  057  845  (9  )  kHz  ,  (1.1) 

6Eth   =  1  057  849  (11)  kHz  ,  (1.2) 

2.  Electron's  anomalous  magnetic  moment 

Oe^p  =  1  159  652  200  (  40)  x  10"^^  ,  (1.3) 

ath  =  1  159  652  460  (127)  x  10~^^  ,  (1.4) 

3.  Fine  structure  constant 

a-ip  =  137.035  993  (5)  ,  (1.5) 

aj^  =  137.035  989  (3)  ,  (1.6) 

where  the  subscripts  "exp"  and  "th"  stand  for  experimental  and  theoretical  values. 
The  quantity  enclosed  in  parentheses  represents  the  uncertainty  in  the  final  digit 
of  numerical  value.  The  value  of  a~^  is  based  on  the  very  accurate  measurements 
of  2e/h  (O.OSppm)  by  the  ac  Josephson  effect.  As  can  be  seen  from  the  results 
given  previously,  theoretical  calculations  using  renormalization  theory  match  the 
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experimental  results  on  the  order  of  O.lppm.   No  other  theory,  to  our  knowledge, 
can  claim  such  a  success. 

If  we  accept  this  remarkable  agreement  between  theory  and  experiment  as  evi- 
dence that  the  fundamental  equations  are  correct,  we  are  led  to  favor  the  conclusion 
that  the  infinities  that  arise  in  the  theory  come  from  an  inadequate  way  of  solving 
the  fundamental  equations.  The  study  of  other  model  field  theories  has  led  to  the 
suspicion  that  the  perturbation  expansions  after  renormalization  become  series  that 
are  only  asymptotically  convergent.  This  has  led  to  a  search  for  a  better  technique 
of  solving  the  equations.  The  technique  should,  in  principle,  allow  the  calculation 
of  the  bare  (noninteracting)  mass,  charge  and  constants  Zi,  Z2  and  Z3. 

In  this  search,  Gell-Mann  and  Low-^^  sought  to  demonstrate  that  the  renor- 
malization constants  are  infinite.  They  stated  that,  although  they  could  not  rule 
out  the  possibiHty  of  infinite  coupling  constants,  its  was  possible  to  isolate  a  nec- 
essary condition  for  which  the  vacuum  polarization  is  finite.  In  a  long  series  of 
papers,  Johnson,  Baker  and  Willey,^^"^^  extended  their  work  and  showed  that  if 
an  eigenvalue  condition  is  satisfied,  then  all  renormalization  constants  in  QED  can 
be  finite.  The  eigenvalue  condition  defines  a  function  that  is  the  coefficient  of  the 
logarithmically  divergent  integral  appearing  in  the  photon  propagator  calculated  in 
massless  QED  due  to  only  those  graphs  with  one  closed  fermion  loop.  This  eigen- 
value condition  is  expressed  in  terms  of  the  bare  couphng  constant.  Adler^^  then 
showed  that  the  zero,  if  it  exists,  must  be  an  essential  singularity  of  this  function 
with  all  its  derivatives  zero  at  the  singularity.  The  existence  of  this  function  has 
not  been  proven.  Nevertheless,  it  is  possible  to  speculate  that  the  existence  of  an 
infinite  order  zero  will  never  be  seen  in  any  finite  order  of  perturbation  expan- 
sion. These  results  generate  great  interest  in  a  number  of  people-'^  for  finding,  in 
a  non-perturbative  manner,  the  solutions  to  the  Schwinger- Dyson  equations.  From 


them,  Zi,   Z2,  and  Z3  could  be  identified  and  it  could  be  determined  whether  or 
not  infinities  are  inherent  in  the  theory. 

Following  this  approach,  H.  S.  Green,  J.  F.  Cartier,  and  A.  A.  Broyles  started 
the  long  term  project  of  solving  the  Schwinger-Dyson  hierarchy  of  equations  using 
a  non-perturbative  method.  The  equations  in  the  hierarchy  are  written  in  terms  of 
momentum  space  variables,  that  is,  in  terms  of  Fourier  transforms.  Figure  1.1  shows 
a  general  block  diagram,  that  can  be  use  in  describing  the  state  of  development  of  the 
project.  Each  block  represents  an  important  step  in  the  solution  of  the  hierarchy. 

Starting  with  the  Schwinger-Dyson  hierarchy  of  equations  (block  1),  H.S. 
Green,  J.F.  Cartier,  and  A. A.  Broyles^° (herein  referred  to  as  Ref.  I)  were  able 
to  determine  the  unrenormalized  electron  propagator  using  an  approximate  form  of 
the  photon  propagator.  In  order  to  obtain  this  solution,  two  approximations  were 
m.ade  (block  2).  Their  first  step  was  to  truncate  the  infinite  hierarchy  with  the  aid 
of  Ward's  identity.^-^  Ward's  identity  relates  the  vertex  function  to  the  next  lower 
equation  in  the  hierarchy  (i.e.  to  the  electron  propagator)  in  the  limit  of  zero  mo- 
mentum transfer  to  the  electromagnetic  field.  The  second  approximation  replaced 
the  photon  propagator  with  its  mass-shell  form. 

With  these  approximations,  they  were  able  to  calculate  the  electron  propaga- 
tor over  the  entire  range  of  the  variables  upon  which  it  depends  (block  3).  They 
found  that  in  order  to  have  a  solution,  the  bare  mass  of  the  electron  must  be  set 
equal  to  zero  and  a  particular  gauge  must  be  chosen,  the  so-called  Landau  gauge. 
These  results  are  in  agreement  with  the  findings  of  Baker  and  Johnson.  ^^  As  an 
additional  result  of  the  determination  of  the  electron  propagator,  the  value  of  the 
renormalization  constants  Zi  and  Z2  were  calculated.  They  conclude  that  within 
the  boundaries  of  the  approximations  made,  Zi  —  Z2  =  1.  No  infinities  appeared 
in  obtaining  the  solutions. 
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Figure  1.1  Block  Diagram  of  iteration  method  for  solving 
the  Schwinger-Dyson  hierarchy  of  equations. 


In  order  to  determine  the  photon  wave-function  renormahzation  constant,  as- 
sociated with  the  photon  propagator  Z3,  it  is  necessary  to  devise  a  more  accurate 
form  for  the  photon  propagator.  The  Schwinger- Dyson  equation  for  it  involves, 
however,  the  electron  propagator  and  vertex  amplitude.  The  electron  propagator 
found  in  Ref.  I  is  available,  but  a  vertex  amplitude  must  be  found.  It  is  possible 
to  demonstrate,  by  using  symmetry  and  invariance  arguments,  that  this  amphtude 
can  be  expressed  in  terms  of  eight  scalar  functions  of  the  electron  and  photon  mo- 
menta. Some  progress  has  already  been  made  on  this  subject.  Using  the  same  type 
of  approximations  which  led  to  the  determination  of  the  electron  propagator,  in  Ref. 
I,  and  a  generalized  form  of  Ward's  identity  derived  by  H.  S.  Green,^^  J.  Cartier  et 
g^^_23,24  (j^gj-ein  referred  to  as  Ref.  II),  determined  all  eight  scalar  functions  needed 
to  define  the  vertex  function  in  the  perturbation  region.  Since  a  direct  analytical 
solution  to  the  vertex  amphtude  is,  at  the  present  moment,  not  feasable  (block  4),  a 
somewhat  different  approach  was  used.  Making  some  reasonable  assumptions  about 
the  behavior  of  the  vertex  function  near  the  mass  shell,  an  approximate  equation  to 
the  vertex  equation  was  derived.  We  call  this  equation  the  model  equation  (block 
5).  The  solutions  to  this  equation  will  indicate  the  functional  form  of  the  vertex 
amphtude  and,  in  principle,  could  be  used  as  starting  solutions  in  a  numerical  it- 
eration procedure.  Also,  a  parametrized  version  of  these  solutions  can  be  use  in  a 
minimization  procedure. 

The  model  equation  was  solved  analytically  by  an  iteration  procedure  and 
tested  to  see  if  the  the  approximations  made  were  valid.  The  solutions  found  proved 
to  represent  reasonably  well  the  vertex  function  near  the  mass  shell  and  to  be  con- 
sistent with  the  Ward  Identity.  Once  again,  no  infinities  were  incurred  in  obtaining 
these  solutions. 

The  next  logical  step  to  follow  is  to  find  the  form  of  the  vertex  function  in 
the  asymptotic  (i.e.    large  momentum)  region  and  which  join  smoothly  with  the 
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perturbation  solutions  found  in  Ref  II  (block  6).  This  is  the  problem  that  we 
address  in  this  work.  The  knowledge  of  the  form  of  the  vertex  function  in  the 
asymptotic  region  will  open  the  possibility  of  determining  the  photon  propagator, 
the  renormalization  constant  Z3  and  thus  the  bare  charge  of  the  electron  (block 
7).  Finally,  to  show  that  the  procedure  is  self-consistent,  at  least  one  iteration  is 
necessary.  With  the  new  vertex  function  and  photon  propagator  a  recalculation  of 
the  electron  propagator  is  possible.  The  present  stage  of  the  project,  including  the 
present  work,  is  delimited  in  the  block  diagram  by  the  broken  line. 

The  present  work  is  arranged  in  the  following  way.  Chapter  II  gives  a  general 
description  of  the  method  used  in  truncating  the  Sch winger- Dyson  hierarchy  of 
equations.  In  this  chapter,  we  also  show  how  to  convert  the  truncated  hierarchy  into 
a  set  of  coupled  differential  equations  and  the  derivation  of  the  model  equation.  In 
Chapter  III,  the  model  equation  is  transformed  into  a  set  of  nine  scalar  differential 
equations.  Also,  a  simplified  version  of  these  differential  equations  is  found  to  be 
valid  in  the  asymptotic  region.  In  Chapter  IV,  the  solutions  to  these  equations  are 
found.  In  the  last  chapter  we  summarize  and  discuss  the  nature  of  the  solutions. 

It  is  felt  that  this  work  can  make  a  positive  contribution  to  a  better  under- 
standing of  the  interaction  between  electrons  and  photons,  and  that  it  will  shed 
some  light  on  the  question:  Is  QED  a  finite  and  mathematically  sound  theory?  In 
addition,  the  techniques  developed  may  prove  to  be  useful  in  the  solution  of  other 
problems  in  QED. 


CHAPTER  II 

DERIVATION  OF  THE  MODEL  EQUATION 

2-1    The  Schwinger-Dyson  Equations 

Quantum  Electrodynamic  Theory  (QED)  deals  with  the  description  of  the 
interaction  of  Hght  and  matter.  Such  interaction  in  its  most  elementary  stage 
studies  the  electron-photon  interaction.  This  description  is  best  handled  mathe- 
matically using  the  concept  of  a  propagator  for  each  of  the  involved  "particles"  and 
an  interaction  site  called  a  vertex.  The  relationships  between  these  propagators 
(Green's  functions)  and  the  vertex  are  expressed  in  a  hierarchy  of  equations  known 
as  the  Schwinger-Dyson  equations.  To  express  formally  these  equations  it  is  nec- 
essary to  define  the  meaning  of  a  propagator  for  each  of  the  "particles"  and  the 
vertex  function. 

The  probability  of  finding  an  electron  at  some  point  in  space-time  given  that  it 
was  at  a  different  point  in  space-time  can  be  computed  from  its  Feynman  propaga- 
tor or  Green's  function.  This  propagator  satisfies  a  differential  equation  analogous 
to  the  wave  function  equation.  In  coordinate  space  the  wave  function  of  a  nonin- 
teracting  electron  ^  satisfies  the  differential  equation 

{if  -?72o)£  =  0  (2.1) 

where  y    =  7'^  d/dx^  and  7'*  are  Dirac  matrices.  The  noninteracting  photon  wave 
equation  is  written  as 

n^A'^   =   0  (2.2) 
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where  D^  is  the  D'Alambertian,  defined  as 


°'    -  ^''-Lt  (2.3) 

and 

r  +1     if  ;U   =   1/   =   0; 

g^""  ^\     0     if/i  7^  z/;  (2.4) 

l-l     if;<   =1/   ^   0  . 

The  electron  propagator  satisfies  an  analogous  equation 

{iiJ   -   m^)Ss,{x',x)    =   8^{x-x').  (2.5) 

Taking  the  Fourier  transform  of  Eq.(2.5)  in  momentum  space  it  is  found  that 
the  noninteracting  electron  propagator  is  given  by 

^^^  =  ^  ■  P.6) 

The  noninteracting  photon  propagator  satisfies  a  similar  equation 

U^DQ{x-x')  =  ib^{x-x'),  (2.7) 

and  the  Fourier  transform  of  the  photon  propagator  yields 

^o(g)  =  -^.  (2.8) 

These  two  solutions  provide  a  complete  description  of  the  electron  and  photon 
when  there  exists  no  interaction  between  them.  When  interaction  is  allowed,  an 
infinite  hierarchy  of  non-homogeneous  integro-differential  coupled  equations  arises. 
The  exact  solutions  of  these  equations  were  not  known.  This  hierarchy  of  equations 
was  formulated  in  their  pioneering  work  by  Schwinger^  and  Dyson'^  and  was  known 
thereafter  as  the  Schwinger-Dyson  hierarchy. 
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Using  the  notation  of  Bjorken  and  DrelP^'^^  these  integral  equations  appear 


as 


or  equivalently 


where 


Sip)  -  5o  +  s{p)  Up)  s{p) 

s-\p)  ^  s^\p)  -  Up) 

Up)  =  j0^  j^\p.  i)S{q)D,^{p  -  q)Y  d\  , 


wik"^) 


le'n 


(2vr: 


Tr 


Y  S{q)r{qA+'k)S{q^k) 


d^q 


r'^(p,3)-7o'  +  A"(P,5), 


and 


L'{PA)  = 


ICn 


(2vr) 


D''\~k')V,ip,p-k)S{p-~k) 


X  T^(p  -k,q-  k)S{q  -  k)TJq  -  k,  q)  d^k 


(2.9) 

(2.10) 

(2.11) 
(2.12) 
(2.13) 
(2.14) 


(2.15) 


+  ■■■  j---  jv\p,p-k)S{p-k).--Y--- 
X  V^{q  -  kn,  kn-i)S{p  -  kn)T^{q  -  k,  q) 
xd^K---d^ki{2T:)-^''  +  ■■■   . 
The  zero  subscript  follows  all  bare  quantities,  that  is,  those  functions  or  con- 
stants which  are  associated  with  noninteracting  particles.    An  overbar  is  used  to 
represent  four-vectors  and  an  underline  to  represent  matrices.    Here  cqA'^  repre- 
sents the  sum  of  all  possible  three-externahpoint,  nodeless  Feynman  diagrams.   A 
node  is  a  bare  vertex  that  if  removed,  leaves  the  diagrams  separated  into  at  least 
two  unconnected  parts.  It  is  possible  to  relate  F''  to  a  four-point  diagram  to  either 
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a  dressed  electron,  or  a  dressed  photon  lines.  These  diagrams  can  be  separated  out 
of  the  four-point  nodeless  diagram  CqE  "  leading  to  the  equation 

x[r^iq,q  +  k)S{q  +  k)r{q  +  k,p  +  k)  +  E^''(q,p-q,p+k)^    d\  .         (2.16) 

The  next  step  in  the  hierarchy  relates  the  four-point  diagram  to  the  five-point 
and  so  on.  In  addition  to  the  Schwinger- Dyson  integral  equations,  there  exist  some 
relations  between  amplitudes  whose  numbers  of  external  points  differ  by  one.  The 
lowest  one  relates  the  vertex  amplitude  V_^  to  the  electron  propagator  5,  and  it  is 
attributed  to  Ward  and  Takahashi.-^^  The  one  relating  the  four-point  amplitude  E^^ 
and  the  three-point  amplitude  F  was  derived  by  H.S.  Green. ^^  Ward's  identity  can 
be  expresses  as  follows: 

{P,  -  g,)r(p,g)  =  S-\p)  -  S-\q)  (2.17) 

or  in  terms  of  A'' (p,  q) 

{p,  -  q.WiP,  q)  =  m)  -  Up)  ■  (2.18) 

An  analogous  equation  can  be  written  for  the  four-point  amplitude  in  terms  of  the 
vertex  amplitude  and  is  given  by 

(pa  +  ^a  -  qx)E^\p.  k.  q)  =  r  (9  ^  k,  q)  -  T (p,p  +  k) 

(2.19) 
^K^{q-k,q)-K'{p,p-Yk) 

and 

k,  E^\p,k,q)  =  T\p-k,q)-T\p,p+k)  .  (2.20) 

Similar  relationships  exist  for  the  remaining  n-point  diagrams.    These  identities 
exactly  define  the  longitudinal  components  of  the  n-point  diagram  in  terms  of  the 
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(n-l)-point  diagrams.  If  we  orient  the  //*''  axis  along  the  n'-^  direction  of  A;,  then  it 
can  be  shown  that  E^"  can  be  written  as 


E^'^iP.  h  q)  =  (P,  +k,-  q,)E^\p,  k,  q){p^  +  k^  -  q^) 


(2.21) 


+  k^E^'ip,~k,q)k-'  +  E^''(p,k,q) 

where  E^  "  is  transverse  to  both  of  the  attached  photon  four- momenta.  Substituting 
equations  (2.19)  and  (2.20)  into  (2.21)  gives 

E^"{p,k,q)  =  \v\p-k,q)-T^{p,p^k)\  (p"  +  F  -  q'') 


+ 


T\p-k,q)-T\p,q  +  k) 


k" 


(2.22) 


lAi/ 


+Er{p,k,q) 


where 


nXi/ 


k,Et''ip,k,q)  =  0 


(2.23) 


Equation  (2.20)  can  also  be  written  in  derivative  form  if  we  orient  the  coordinate 

axis  so  that  the  i/*''  axis  lies  along  k.   Then  dividing  by  k"  and  taking  the  Umit  as 

k"  vanishes  gives 

dT\p,q)        dT\p,q) 

dpu 


i\u 


E^'iPAq) 


dqu 


(2.24) 


Similarly,  letting  kx  approach  {qx  -  px)  in  equation  (2.19)  gives 

dV\p,q)       dT\p,q) 


-[Xu 


E  "{p.q-p.q) 


(2.25) 


dpv  dqi, 

Using  similar  arguments,  it  is  possible  to  separate  V"  into  longitudinal  and  trans- 
verse parts  so  that 


L'iP.  q)  -  [Sr\p)  -  S-\q)]    k^  +  rf  (p,  q) 


(2.26) 


where  k^T!^  =  0  or  explicitly 


r.iq  +  k.p+k)  -  [Sr\q  +  k)-S-\p-{-k)] 


{q'-pn 


(q~p) 


fi^2 


(2.27) 
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Using  equations  (2.16),  (2.21)  and  (2.27)  leads  to  the  following  equation  for  the 
vertex  amplitude: 

T\p,p+  k)=l^  +  J^   f  D,,{p  -  qh''Siq)\T\q,  q  +  k)S{q  +  ~k) 


(27r)^ 


r;{q  +  k,p  +  k)  +  [s-\p)-s-\q)] 

[r-^^q  +  k,p+k)~r{q,p)]    k^ 


{q  -  P? 


(2.28) 


+ 


V\2q~p,p^k)-T\q,2p^k  -  q)     (p"  -  q'') 


+  Ef(q,p-q,p+k)jd'q. 
2-2    Approximations  to  the  Schwinger-Dyson  Equations 

Although  equation  (2.28)  is  exact,  it  involves  another  equation  in  the  hierar- 
chy through  the  four-point  amplitude  E^" .  It  is  possible  to  close  the  hierarchy  of 
equations  by  noticing  in  equation  (2.28)  that  the  integrand  has  a  maximum  at  the 
pole  {p  —  q)  —  0  and  therefore  it  is  reasonable  to  approximate  E^ "  by 


Xv, 


Er{q,p-q,P  +  q)  -  Ef{q,0,p  +  q)  . 


(2.29) 


Since  A  is  transverse  to  A;,  and  u  is  transverse  to  (p  —  q),  E_t^  can  be  written  as 

^^  _  k^^r^u  _  {Pv  -  (JnXP"  -  ^'') ^Xr,  _ 


771AJ/  rpXi 


k^ 


ip-<l) 


2  — 


(2.30) 


Substituting  equations  (2.24)  and  (2.29)  into  (2.30)  gives 


E^''{q,p-q,p  +  q)-  -[g^g'a 


k^k") 


.  _  APc.-qa)ip''-qn.dr>iq,q  +  k) 


{p  -  qf 


dqa 


(2.31) 


which  can  be  used  to  eliminate  E^"  in  equation  (2.28).  This  equation  coupled  with 
equations  (2.9)  to  (2.16)  could  be  solved  in  principle  for  F   ,5  and  D^^. 
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It  is  possible  at  this  stage  to  solve  equation  (2.28)  for  the  component  of 
E  iPiP  +  k)  transverse  to  k^  and  then  use  equation  (2.26)  to  construct  T^.  Tak- 
ing the  transverse  components  of  F^  will  ehminate  all  terms  proportional  to  k-^. 
Substituting  equation  (2.31)  into  (2.28)  gives 

T^iP,  p+-k)  =  j^^  +  ^jD,,{p-  q)2''Siq){L^iq,  q  +  ~k)Siq  +  k) 


X 


+ 


T-;iq  +  k,p+-k)+[S-\p)-S-\q)]     ^-J^^ 


Tt(2q-p,p  +  k)  -  Tf{q,2p+k  -  q)     {p^  -  q'') 


I    drXq,q  +  k)                                        drf{q,q  +  k)]   ,, 
+  g^^ ip.-q.)ip    -q) ^^- 1^5- 

(2.32) 
The  final  goal  would  be  to  solve  this  equation  with  Eqs.(2.9)  to  (2.13)  self- 
consistently.  To  this  end,  it  will  be  advantageous  to  solve  a  much  simpler  equation 
obtained  by  making  some  seemingly  sound  approximation.  The  solutions  to  this 
equations  will  shed  light  on  the  form  and  behavior  of  the  complete  solution  in  both 
the  perturbation  and  asymptotic  limit. 

The  last  step  in  the  procedure  of  reducing  the  Schwinger- Dyson  equations  to 
a  more  tractable  form  was  to  convert  the  integral  equations  into  a  set  of  differential 
equations  with  appropriate  boundary  conditions.  This  method  was  first  developed 
by  H.  S.  Green  in  connection  with  the  Beth-Salpeter  equation.  It  was  first  used  in 
the  study  of  the  Schwinger-Dyson  equations  by  Bose  and  Biswas. ■^^ 

2-3   The  Model  Equation 

It  is  possible  to  recast  equation  (2.16)  using  Ward  and  Green's  identities  for 
both  terms  in  the  square  brackets.  Replacing  p  hy  q  and  making  use  of  Ward's 
identity  gives 

r\g-  q  +  k)Siq  +  fc)r,(g-  +  fc, p  +  fc)  ^  V\q,  q  +  k)S{q  +  k)  ^~    '^^'^  ^^   .     (2.-33) 
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Substituting  this  equation  and  the  expression  for  E^^{q,  0,q  +  k)  in  equation  (2.16) 
and  taJjing  the  transverse  component  yields 

TXp,P+  k)  =  j^^  +  -^J  D.^ip  -  q)Y^F'''{q,  q  +  k)d\  (2.34) 


where 


F^%q,q  +  k)  =  S{q)~  {r^{q,q  +  k)S{q  +  k)}  S'^q  +  ~k)  .  (2.35) 

Using  the  renormahzation  group, ^^  arguments  have  been  presented  to  show 
that  the  asymptotic  form  of  D^^  is  the  same  as  that  near  the  mass  shell.  One  of  these 
arguments  is  that  the  contribution  of  the  vacuum  polarization  to  the  Lamb  Shift  is 
only  27  megahertz  out  of  1058.  The  vacuum  polarization  gives  a  measurement  of 
the  departure  of  the  photon  propagator  from  its  mass  shell  form.  It  is  also  observed, 
once  again,  that  the  integrand  of  equation  (2.34)  is  largest  when  the  argument  of 
D^i,  vanishes.  Hence,  the  photon  propagator  is  approximated  by 


D^M 


9fiu    1        _2 


(2.36) 


Q     J  1 
where  Z3  is  the  photon  renormalization  constant. 

The  gauge  is  chosen  to  be  that  found  necessary  to  obtain  a  finite  solution  to 
the  electron  propagator  equation^"  with  a  vanishing  bare  mass  for  the  electron.  In 
the  above  reference  the  electron  propagator  is  found  to  be  given  by 


S-\p)^^Aip')+B(p')P 


where 


A{p^)  =  m 


3a(rn   —p   ^/Ai^p 


B{p^)  ^  1 


(2.37) 

(2.38) 
(2.39) 
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and  a  ~  1/137  is  the  fine  structure  constant.  With  these  approximations  it 
is  possible  to  transform  equation  (2.32)  to  a  differential  equation  by  taking  the 
D'Alambertian  in  momentum  space  (see  Eq.(2.3))  on  both  sides  of  Eq.(2.34).  Us- 
ing the  identity 


d 


it  is  possible  to  write  Eq.(2.36)  as 


D,.{p-q)  =  -\z,^^ 


HP  -  qf  =  2 


Siiu 


(P  ~  qf 


^ip, -q,)ip.- q.)  (2_4Q) 

{p-qy  ^ 


'Zs-       ' 


d 


2       (p  -  qf       4       d{p^  -  qf,)  d(p^  -  q^) 
Application  of  the  D'Alambertian  to  Eq.(2.41)  gives 


\n{p-qY  (2.41) 


U^D,.{p-q)  =  -\z,g,,U'j^^  -  \z,d,dj[}^\n{p-qf 


Substituting  the  identities 


(2.42) 


D'(--)I-^'(2-)'^'(^-^-)' 


n^  ln(p-  -  qf  = 


{p  -  qf 


into  Eq.(2.42)  gives 

^'D,,{P  -  q) 

Therefore, 


-z. 


27r^  i  QiJiuS'^ip  -  q)Zz 

d  d 


Q(^pii  _  qti)  Q(^pu  _  qu-^  (^  _  ^y 


n'Ttip,p  +  k)  = 


te 


(Stt^) 


ie^.Z: 


-2iTHZ,-fS{p) 


d 
dPn 


(2.43) 
(2.44) 


(2.45) 


(27r) 


0^3      ,    1^,, 


Lt(P,P  +  q)S{P  +  q)    S~\p  +  q) 
d  1 


d{pf-  -  Q^)  dip"  -q''){p-  qf 


d 


xS(q)^    r^{q,q  +  k)Siq  +  k)    S-\q  +  k)\d'q 


(2.46) 
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In  terms  of  the  tensor  F"^" ,  Eq.(2.43)  becomes 


n'rXp,p  +  k)  =  lejF'%p,p  +  k)-  ''     ^      ^ 


2  -^- 


/ 


1  ^A. 


(27r)4  J    d{p^-qu)      {p-qf 


ll  F^'  +  G^ 

2-Lv 


where 


G^  = 


d 


J  -  d{p-  - 


(27r)2  dp"  J  -    d{p°  -  q") 
Applying  ^  to  G^  and  using  Eq.(2.43)  we  obtain 


dp" 


and  hence 


(P-I) 


F^'d^q 


(2.47) 


(2.48) 


(2.49) 


y^n^T\p,p  +  k)  =  e 


lfl^F'-'{p,P  +  k)  +  ^F''^{p,p  +  k) 


(2.50) 


We  have  define  e^  =  egZa  and  e  =  e^/(27r)^  =  {oi/^)  where  a  again  is  the  fine- 
structure  constant.  It  is  possible  to  write  Eq.(2.50)  in  a  more  symmetrical  form 
as 

□'r\pi,P2)  =  e 


\T^F^'{pi.P2)+^y-'F^\p,.P2) 


where 


and 


:.\u 


A/-     - 


£-^"(^1,^2)  =  -S{p^)V\p^,p,)Sip2)r  +£\puP2) 


(2.51) 


(2.52) 


/A/- 


A"    (Pl,P2)  =  S{pi) 


dT\puP2: 
dpu 


r\p„p2)S{p2: 


dS-\p2)        A 
dp.  ^  , 


(2.53) 


Here  pi  and  p2  refer  to  the  outgoing  and  incoming  electron  momenta  respec- 
tively. The  tensor  F  "  has  been  written  in  such  a  way  that  it  supports  further 
approximations.   To  first  order  in  e  it  is  well  known  that  near  p^  =  p^  =  rr?  (i.e. 
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near  the  mass  shell)  T^  ^  7^ ^2^  ^^^^  therefore  ^^^P^  =  X^^2^-  ^^  ^^  as- 
sume that  this  form  is  correct  we  can  neglect  the  derivatives  of  the  vertex  function; 
^•e.  -gpr  ^  0.  Notice  also  that  to  this  order  of  approximation  Zi  ^  Z2  !^  1.^°  Under 
these  approximations  it  is  reasonable  to  assume  that  A^"  =  0.  This  leads  us  to  a 
simpler,  more  compact  equation  which  we  called  the  model  equation  for  T^  and  is 
explicitly  written  as 


n^T\p„p,)  =  -e 


Ix^  A\puP2)  x"  +  f-'A^f 


where 


A\puP2)  =  Sip^)T\p^,p2)S(p2: 


(2.54) 


(2.55) 


We  are  interested  in  the  solutions  to  this  equation  coupled  with  equations 
(2.36)  and  (2.37)  in  the  asymptotic  region.  These  solutions  must  join  smoothly  with 
solutions  near  the  mass  shell  (i.e.  perturbation  hmit).  To  solve  these  equations,  it 
is  convenient  to  write  them  in  term  of  invariant  functions.  This  can  be  done  as 
follows:  let  us  define  the  matrices  (see  Appendix  A) 


5^,P 


1,.  -  [l^'lj         1^^^  =  I  hyl.J  =  ^ex.^px'r 


(2.56) 


I5 


-^7o7i72l3 


(2.57) 


where  £0123  =  -£°^^^  =  1  .  Here  the  symbols  [,]  and  {, }  represent  the  commutator 
and  anti-commutator  of  Dirac's  matrices.  One  can  also  define  the  following  vectors 
and  tensors 


C^  = 


r^  - 

^liu  — 


\Tr 

pA 

5 

\Tr 

1 

\1' 

r  7 

\lr 

VK 

-> 

(2.58) 
(2.59) 
(2.60) 
(2.61) 


The  tensor  equations  derived  from  Eq.(2.54)  are 
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D^^^Tr 

A^' 

1 

(2.62) 

fi          4 

^^'l,J 

) 

(2.63) 

-^/ii/        4 

^^'l,' 

^     , 

(2.64) 

D^      =  ^  Tr 

A^; 

^.^ 

,p      ~  ^^//f  per-'  '^ 

^\a' 

(2.65) 

u^c^  =  _3er>' 


D^C^  =  2eDj  -  2ea^D-2a''D^  , 


-tDl,  -  led^U-'^d'Dlp  +  2ea,n-2c/''£»i,  , 


-26a^n-2a'^i),\, 


26a,n-2a-Dj^, 


2e5,D-2a-Dj,^ 


(2.66) 
(2.67) 
(2.68) 


(2.69) 


If  we  write 


Skfx 


^1 


Sin)  = 


(p1  -  A\)  ' 

where  Ai  =  Ai(p^),  A2  =  ^2(^1),  and  also  define 


{Pl-AD 


(2.70) 


D^{pl-Al){pl-Al), 


(2.71) 


then  equations  (2.58)  to  (2.61)  become 

D^  -  \Tr[T\p,,p^){p2^A2){p,^A,)\/D 

=        [{A,A2  +  pi  •  P2)C^  +  (A1P2'  +  A^p>t)Cl 
-{ptP2)CUlD, 
Dl  =  \Tr[v\p^,p,){h^A^h{^^-vA,)\lD 


(2.72) 
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[{A2Pl^,+AiP2^)C^+(AU2-pi  ■P2)C^ 

+  {Pi,p'2  +  P2MCt  +  {A^pl  +  A,p'^)Cl, 
-{PiP'2K.,]/D.  (2.73) 

iPj,  =  iTr[T\puP2)ip2+A,)j^^if>,  +  A,)]/D 

=  [--iPliiP2u  -  P2iiPlv)C^  +  {A2Pxu  -  AyP2u)C^ 

-  {A2Plt,  -  AiP2^)Ct  +  {A1A2  +  Pl  ■  P2)C^. 

-  {PlvP'2  +  P2.p'i)C^p  +  iPl^p'2  +  P2MCtp 

+  (A2P?  +  ^P20CiJ/^,  (2.74) 

£>^p  =  \Tr[T\p,,p2){p2  +  A2)x^^^{p,^A,)\lD 

=  [(Pl/xP2t/  -  P2i,Plu)Cp   +  ipiuP2p  -  P2vP\p)Cl 

+  {PlpP2p.  -  P2pP\^)C^  +  {A2P\p  +  AiP2m)C^p 

+  (A2Pi^  +  A^p2^)C^p^  +  (A2P1P  +  A^p2p)Cl^ 

+  (P1.P2  +P2.PnC^p\.  +  {PipP2+P2pPl)Ci^.MD  .  (2.75) 

In  deriving  the  above  expressions,  several  matrix  identities  were  used  (see 
Appendix  A).  It  is  possible  to  express  these  tensors  in  terms  of  scalar  functions 
using  the  following  reasoning.  To  construct  the  vector  fiinction  C'^  there  is  only 
one  vector  at  our  disposal,  p  since  the  vector  k'^  does  not  appear  in  the  transverse 
part  of  r  .  Therefore,  the  most  general  form  of  C^  is 

C^  =  up^  .  (2.76) 

To  construct  the  tensor  C^  we  have  at  our  disposal  only  the  tensors,  8^,  P\p.p^-,  and 
P2p.P^-  Hence,  the  most  general  form  for  C^  is 

Ci  =  v8^^  +  {u2Pi^  +  Uip2^)p^  .  (2.77) 
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In  constructing  the  general  form  of  C^^,,  we  notice  from  Eq.(2.60)  and  the  antisym- 
metry property  of  7  ^  under  the  exchange  of  the  indices  that  C^^  —  -C^  .  The 
only  tensors  at  out  disposal  to  construct  C^^  are:  pi^S^,  PiJ^,  P2^^^  P2Ji  and 
PiiJiP2^P^-  Thus,  the  most  general  form  of  C^^  which  satisfies  the  antisymmetry 
relation  is 


C^u  =  {V2PI,^  +  ViP2u)Sl  -  iv2Plf,  +  ViP2f,)~6^ 

+  U3{pi^tP2u  -  P2^Plu)P^  ■  (2.78) 


Similar  arguments  can  be  used  to  conclude  that 

C^up  =  -"Z  [{PlvP2p  -  P2uPlp)~S^  +  {pipP2^  -  P2pPl^i)~^t 

+  {pip^P2u  -  P2p.Pi^)~S^]  .  (2.79) 

With  these  definitions  the  vertex  amplitude  takes  on  the  form 

£  {Pi,P2)  =  1_^v{p\,pI,p^-P2) 

+  p^u(pi,P2,Pi  -Pi) 
+  ^^P^Ui{pI,pI,pi  •p2) 
+  ^2P^U2{pI,pI,Pi  ■P2) 
+  P^[^i^^2\'^i{p\^pl^Pi  ■P2) 
+  [/i,7^]«i(p?,pIpi-P2) 

+  [^2'7^]"2(Pi,P2,Pl  ■P2) 

+  ie^"'%PiuP2pV3ipl,plpi  ■P2)  .  (2.80) 

A  tilde  is  used  over  variables  to  specify  transverse  parts  to  k.  Exphcitly,  the  trans- 
verse part  of  p^  is  given  by  p^  =  p^  -  ^k^.   Also  the  tensors  D^,  D^,  D),„  and 

ft  til  '  LLl' 
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D^„^  can  be  written  in  terms  of  scalar  functions  as  follows: 


Di-=sS^  +  (R2Pl^  +  RiP2^)p^  , 
K  =  (P2  j;  -  P2,)S^)S,  +  {p,J^,  -  Pl,)~6^)S2 

+  R3{pi,iP2u  -  P2,iPlu)p      , 
D^up  =  •^a  [iPluP2p  -  P2vP\p)Sl  +  {pipP2n  -  P2pPlti)^^ 
+  {PlnP2u  -P2^lPlu)Si]   ■ 


(2.81) 
(2.82) 

(2.83) 

(2.84) 


Substituting  these  expressions  into  equations  (2.72)  to  (2.75)  gives  the  following 
results  for  the  scalar  functions  R,Ri,. . .  etc.: 


R 


Ri^ 


R2  = 


R, 


5  = 


Si 


(A1A2  +Pl  •  p2)u  +  {A1P2  +  A2P1)  ■{uiP2  +  U2P1) 

+  (Ai  +  ^2)^  +  (pI  -  pi  ■  P2)v2  -  (pI  -  Pi  •  P2)vi 

-\PIPI-(P1-P2f]u3\/D   , 
Aiu  +  A1A2U1  +  pIu2  +  U  +  (Ai  -  A2)ui  -  (A2P1  +  AiPi  ■  P2)uz 


(2.85) 


-{p\-pi-p2)v,^lD,  (2.86) 

A2U+P2"1  +  A1A2U2  +U  +  (Ai   -^2)^2  -{AipI  +  A2P1  ■P2)UZ 

+  (P2  -Pi  -^2)^3  \lD  , 

-U  -  A1U2  -  A2U1  +{AiA2  +Pl  •P2)uz 
■^  {V2  -  Vi)  +  {Ai  +A2)vA/D  , 

(A1A2  +Pl  ■P2)V  +  {A2p\  +  Aipi  ■P2)V2  -  AipI  +A2P1  ■P2)vi 


(2.87) 


(2.88) 


+  {pipi-{Pi-P2r)us\/D  , 


-Aiv  +  A1A2V1  -  p\v2  -  (A2P1  +  AlPi  •  P2)vz  \ID 


(2.89) 
(2.90) 


26 


52  =1  ^2^  +  A1A2V2  -  plvi  +  {A^pI  +  A2P1  ■  P2)V3  \/D  ,  (2.91) 


^3  =<^  U  -  A2V1  +  A1U2  +  (Ai A2  +  Pl  ■  P2)V3  \/D 


(2.92) 


CHAPTER  III 

DERIVATION  OF  THE  SCALAR  DIFFERENTIAL  EQUATIONS 

3-1    Definition  of  Variables 

To  obtain  the  scalar  differential  equations  restilting  from  the  model  equation, 
let  us  define  the  following  variables  which  will  prove  to  be  very  useful  in  describing 
the  asymptotic  region.  Let 

pjPl  ;  (3.1) 

^=U-Tf;         y  =  -iz  +  z-')-  (3.2) 

P2  ^ 


z  =  y  +  Vy^  -I;        z  1  =  y -  ^y2  _  1  _  (33) 

The  inverse  transformations  are 


p{=f  ^2p-h^k'  =xz;         pI=p^ -2p-k  +  k'^  ^xz-^  ,  (3.4) 

so  that 


p^  ^xy-P;      p-k  =  ^x^/y^-l;      pi  ■  P2  ^  p^  -  P  =  xy  -  2P  .  (3.5) 

If  both  pj  andpl  are  time-Hke  or  space-like  4- vectors,  the  variables  x  and  y  are  real. 
With  these  restrictions,  the  variables  y  or  z  are  always  greater  than  one.  Using  this 
choice  of  variables,  the  momentum  space  is  divided  into  three  regions: 

perturbation  region,  x  ?»  m?  ,    y  Ri  1,  P  ?a  0  ; 
inner  asymptotic  region,  x/P  >  m^  ,    y  ~  1  ;  (3.6) 

outer  asymptotic  region,  x/P  :§>  m^  ,    y  >  1  . 
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In  terms  of  these  variables  the  differential  operator  5^  =  djd-p^  can  be  written  as 

a,.  =  i,|:  +  M,|.  (3.7) 

where 

^A.  =  2  {y  p^  -  Vy'-l  ^v}   ,  (3.8) 

M,  =  -2^^^^{v'y2-l  p^-yA;^}   .  (3.9) 

If  we  have  a  vector  function  W^'^  of  the  form  W^  =  w(x,  y)p'^.,  then 


d,w^  =  id,w)p^  +  wl"; 


(3.10) 


and 


U'^W^=p^U^w  +  2d^w  ,  (3.11) 

where  D^   =  d'^d^  is  the  D'Alambertian  operator.     To  specify  the  derivatives  of 
w{x,y)  with  respect  to  the  x  and  y  variables  we  use  the  notation 


dw 
dw 


t^x   -    ^^  ,        W,,  -    ^^^  , 


d^ 


w 


^y  ~  a   '     ^yy 


dy  '      ^^        ay2  ' 


Wxy  — 


dxdy 


(3.12) 


Using  these  symbols,  it  is  found  that  the  D'Alambertian  operator  acting  on  a  scalar 
function  of  the  variables  x  and  y  gives 


U^w{x,  y)  =  4(a;y  -  k'^)wxx  +  4(2y  -  — )u;^  + ^^- '-w 

X  x^ 

P        (y^  -  1) 


yy 


+  4 


y 


Wy  -  4(y    -  l)w 


xy 


(3.13) 


With  the  help  of  Eq.(3.10)  it  is  easy  to  see  that 


W^W^ix,  y)  =  4(,Ty  -  k'^)iVx^  +  4(3y )iu^  + 


4P(y2-l) 


w 


x^ 


yy 


+  4 


■     P       2(y^  -  i: 

^        2 


ty 


J,  -  4(y2  -  l)'u;^3, 


(3.14) 
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when  acting  upon  a  vector  function,  since  from  Eqs.(3.7),  (3.8)  and  (3.9)  gives 


d^w  =  2 


jy'  -  1) 

y  w^ Wy 


P 


(3.15) 


To  distinguish  between  the  two  differential  operators,  we  used  the  subscript  s  or  v 
to  indicate  the  operator  used  in  Eq.(3.13)  or  (3.14),  respectively.  Explicitly, 
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k\  d 


Ul  =  A{xy  -  k")—;  +  ^{2y  - -) 


+ 


5x^ 
4P(y^  -  1)  a^ 


+  4 


X    dx 

P       2(y2  _  1) 


x^ 


X 


d_ 
dy 


-4(y^-l) 


dxdy 


(3.16) 


and 


D^.  =  4(^^-^^)^  +  4(3y-^)^ 


+ 


dx"^ 

4P(y^  -  1)  d^ 
x'^  dy'^ 


+  4 


■    k^      (y'  - 1) 

y 

x^ 


X       \  dy 


d 


-4(y^-l) 


d^ 


dxdy 


(3.17) 


To  recast  the  nine  coupled  differential  equations  that  arise  from  the  model 
equation  in  terms  of  the  new  variables,  let  us  define  the  operators  h  and  /  by 

d     (y2  - 1)  a  \ 


h=   2     y 


dx 


l  =  -2^y2-l 


X        dy 

d_  _  y_d_ 
dx      X  dy 


and  also 


4  =  A  ±  vZZli. 

dx  X        dy 


(3.18) 
(3.19) 

(3.20) 


Notice  that  with  this  notation 


d^  =  p^h 


(3.21) 
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It  is  easy  to  show  that  in  terms  of  these  differential  operators 

dn  =  Pucz-^U  +  P2,izf-  =p^,h^  kj  (3.22) 


since 

z-'u-zf_  =  i. 

3-2    Derivation  of  Differential  Equations 


(3.23) 


At  this  moment  we  turn  to  the  problem  of  transforming  Eqs.(2.66)  to  (2.69) 
into  scalar  differential  equations.  It  is  easy  to  see  that  Eq.  (2.66)  can  be  written  in 
terms  of  the  scalar  functions  u  and  R,  as 

Ulu  =  -3eR  .  (3.24) 

To  obtain  the  scalar  differential  equations  from  Eq.(2.67)  it  is  observed  that  this 
equation  is  equivalent  to  two  coupled  differential  equations.  To  this  end  we  define 
the  vector  function  E^  =  ep^.  With  this  definition,  equation  (2.67)  becomes 

n'Ci  =  2e[D^^-d,E^]  (3.25) 

coupled  with 

D^^;^  ..:  d^D^  .  (3.26) 

It  is  important  at  this  moment  to  reahze  that  from  Eq.(2.67) 

a^C^  =  0  .  (3.27) 

This  follows  from  a  direct  application  of  S''  to  C^.  The  proof  follows: 

=  2e  [d^D^^  -  d'^D^]  =  0  . 
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Assuming  that  d^C^  is  finite  at  p^  =  0,  and  p^  =  oo,  then  the  only  choice  possible 
is  the  result  stated  in  Eq.(3.27).  In  addition  to  this  argument,  we  find  that  the 
tensor  C^  is  defined  in  terms  of  the  scalar  function  v.  In  the  asymptotic  region,  the 
largest  contribution  to  C^  comes  from  this  function.  It  is  found  (see  appendix  D) 
that  in  the  asymptotic  region  the  function  v  approaches  a  constant,  thus  supporting 
the  hypothesis  that  9^C^  ~  0  to  leading  order  in  the  variable  x. 

If  we  apply  in  general,  d'^  to  C^  we  obtain 


(3.28) 


Evaluating  Eq.(3.28)  at  a  =  /i  gives 


d^C^  =  d\  +  5(U1  +  U2)p^  +  p^  [pi  ■  UU2  +  P2  ■  Uui]  =  0 


(3.29) 


or  using  Eqs.(3.21) 


h(v)  +  5(ui  +  U2)  +  pi  ■  Uu2  +  P2  ■  Dui  =  0  . 


(3.30) 


Applying  da  to  Eq.(3.28)  gives 


D^C^  =  ~6^  piv  +  2{u,  +  ti2)]  +  2  [pi^a\2  +  P2^5^"l] 

+  [pi^n^U2  +  P2^^lui]  p^  +  2(a^Ui  +  d^U2)p^ 


(3.31) 


Using  Eq.(3.21)  we  obtain 


U'C^,=6l[Ulv  +  2{u,^-U2)\ 


+  Pi^P 


+  P2^P^ 


Ulu2  +  2h{u2)  +  2z-^U{ui  +  U2) 
Ului-V2h{ui)+      2zf-{ui+U2) 


(3.32) 
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U'Cl=~8l 


Ulv-2h{v)-^2{ui+U2) 


+  P2^f 


nlu2  +  2z-^U{u^+U2) 
Ului+         2zf_{Ui+U2) 


(3.33) 


From  Eq.(2.82),  it  follows  that  the  right-hand  side  of  Eq.(3.25)  becomes 


2e  [Dl  -  d.E'']  =  2eUs  -  e)S^^  +  p^,p^[R2  -  z'' f+i^)] 

+  p,^p^[R,-      ^/_(e)]|.  (3.34) 


Also,  from  a  direct  application  of  5''  to  D^  we  obtain  that 

d'D^  =  [5{R,  +  R2)  +  K5)  +pi  -UR^  +P2  ■  UR2]   . 


(3.35) 


Comparing  Eqs.(3.33)  and  (3.34),  it  becomes  evident  that  we  can  extract  four  dif- 
ferential equations  for  the  functions  u,wi,U2  and  e.  They  are 


Ulv  -  2h{v)  +  2{uy  +  U2)  -  2e(5  -  e)  , 
Ulu,  +      2z/_(ui  +  U2)  =  2e[R,  -      zf^ie)]  , 
Ulu2  +  2z-^U{ui  +  U2)  =  2e[R2  -  z'' U{e)\  , 


and 


Ule  =  5(i?i  +  R2)  +  Pi  ■^R2  +  P2  ■  ^Ri  +  h{S) 


(3.36) 
(3.37) 
(3.38) 

(3.39) 


To  obtain  the  differential  equations  for  the  functions  V\^V2  and  M3,  we  must 
first  simplify  Eq.(2.68).  This  can  be  accomplished  if  we  realize  that 


U'd^Cl  =  ed^D^^  , 


(3.40) 


which  follows  from  a  direct  application  of  d"  to  Eq.  (2.68).  This  implies  that 


eU-^d'^Dl^^d'^Cl, 


(3.41) 
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and  therefore  Eq.(2.68)  becomes 


□'C/,t  +  2{d,d''Ct,  -  d^d^Ci^)  =^  -eDl^  .  (3.42) 

Using  the  definition  of  C^^  given  in  Eq.(2.78),  we  can  write  C^^  as 

Cop  =  {PlpV2  -^P2pVl)~8l  -  {piaV2  ^-p2aVi)8^ 

+U3{piaP2p-P2aPlp)p^   ■  (3.43) 

Applying  d^  to  Eq.(3.42)  gives 

d'C^p  =  [3(t;i  +  ^;2)  +Pl  •□^2  +P2  -a^l]^^  -  {Plad\2  +  P2ad\,) 

+  [PlaP2  ■  ^U3  -  p2aPl  "  ^^3]  p^  +  4:U3{pia  -  P2a)p^   ■  (3.44) 

In  general,  if  we  apply  d""  to  C^^,  this  results  in  the  identity 
d'^C;:,  =     <5?(^i  +i;2)  +  (pi.a%2+P2.a%i)^^ 

+  {PlliP2v  -p2iiPlv)p^d'^Ui  -^Uz{pi^P2u  -  P2tiPlu)S^°' 

+  u,  [^>,  +  Sy,^  -  S^py,  +  <5>  J  p^  .  (3.45) 

Applying  da  to  Eq.(3.45)  and  factoring  all  scalar  functions  multiplying  the  tensors 
iPiJil-Piii^i)^iP2Ji  -P2ii^),  and  {pif,p2v  -P2nP\u)p^  gives 

^''Ci,  =(Pi.^i  -  Pi^8l){Ulv._  +  2z-'U{v,  +  v^)  -  2^3} 
+(P2.^,!  -  P2j^){Ulvy  +  2z/_(t;i  +  ^2)  +  2U3] 
+iPii^P2u  -  P2p.Pi^)p^{Ulu3  +  6h(u3)}  .  (3.46) 
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An  analogous,  but  otherwise  tedious  calculation  yields 

+  (P2Ji  -  P2ii6u)G2{vi,V2,  Us) 
+P^{.PliiP2v  -  P2iJiPlu)Gz{Vi,V2,Uz)   ,  (3.47) 

where 

Gi{vi,V2,us)  =\-^z-^  f^{vi  +V2)-2h{v2) 

^2(^1,^2,^3)  ^l  -     3zf-{vi  +V2)  -2h{vi) 

-      2r/_[pi  •DU2]  -      2/_[p2 -n^^i] -Pi -nua -4u3  L 
^3(^1,^2,^3)  =j  \h[h{v2  -vi)]  -p  -Uuz  -7h{u3)  i  .  (3.48) 

Substituting  Eqs.(3.46)  and  (3.47)  into  Eq.(3.42)  gives 

U'Ci,  +  2{d^d^Ct,  -  d.dfCl^)  =  (p,J^^  -Pi,Si)H,{v^,V2,u,) 

+  {P2uS^  -p2,iK)H2ivi,V2,U3) 
+P^{PltiP2u  -  P2^Pu)H3{vi,V2,U3)   ,        (3.49) 

where 

Hi{vi,V2,U3)  ^<Ulv2  -4z~V+(^i  +^^2)  -4/1(^2) 

-  2z-'f+\p^  ■  nt;2]  -  2z-^f+\p2  ■  Uvi]  +  2p2  ■  Dus  +  6uA  , 
H2(vi,V2,U3)  =lQlvi  -      4:zf-(vi  +V2)  -4:h{vi) 

-  2z/_[pi  •  0^2]  -      2zf_\p2  •  Dui]  -  2pi  ■  n«3  -6uA  , 
H3(vi,V2,U3)  =<Ulus  -  lOhim)  -  2p  •Du3  +  h[h(v2  -  ^i )]  I  .  (3.50) 
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Comparing  Eqs.(3.42),  (3.50)  with  Eq.(2.83)  and  factoring  all  terms  proportional  to 
the  tensors  {piJ^^  -  Pi^~6^),  {p2u~^i  -  P2^^^)  and  p\pn,P2u  -  P2^Pi.)  yields  the 
desired  scalar  differential  equations  for  the  functions  t;i,U2  and  U3.  They  are 

uivr  -  4/z(wi)  -    2zf-\pi  ■  0^2]  -    22/- [p2  -n^i] 

-  2z/_(t;i +U2)-6u3-2pi -D^xa  =2e5i  ,  (3.51) 

n2„u2  -  4%2)  -  2z-7+[pi  •  n^'2]  -  2z-'f+\p2  ■  a^i] 

-  2z-V+("i  +  ^2)  +  6u3  +  2p2  •  Dwa  =  2e  52  ,  (3.52) 
nlus-10h{u3)-2p-Uh{m)  +  h[h{v2-vi)]  =-eR3  ,  (3.53) 

where  the  functions  5i,   ^2  and  R^  in  the  right-hand  side  of  Eqs.   (3.51)  to  (3.53) 
are  given  by  Eqs.(2.88),  (2.90)  and  (2.91). 

To  obtain  the  final  differential  equation  for  the  function  U3,  we  first  simplify 
Eq.(2.69)  by  reaHzing  that 

^'d'C'^up  =  -26^"^;,,  .  (3.54) 

This  equation  still  contains  all  the  information  of  the  original  equation.  It  is  possible 
to  reduce  this  equation  to  a  vector  equation  by  multiplying  Eq.(3.54)  with  the  tensor 
pj'p^.  Due  to  the  antisymmetry  of  C^,^  and  D^.p  under  the  interchange  of  two  of 
its  indices,  and  the  fact  that  C^^^  =  0  if  any  two  indices  are  equal,  we  can  solve  the 
equivalent  equation 

□VrP2a-C,\,  =  n^d"  [p',p;Ci.,]  =  -2ed''  [p^P^,DiJ   .  (3.55) 

Using  Eq.(2.79)  and  contracting  C^^p  with  p'^p^  gives 

PlP2Clup  =  ^3  I  (Pl  •  p2P2p  -  pIpXp)P^  +  (Pl  •  P2PIP  -  P\P2p)P^ 
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+\pIpI  -  iPi  ■  P2)']Si'^  .  (3.56) 

Applying  d^  to  Eq.(3.56)  gives 

d'  [p'iP2Ci.,]  =  -SPvsp^  .  (3.57) 

Similarly,  from  Eq.(2.84)  it  follows  that 

PlPl^lup  ^  SA  (pi  •  p2P2p  -  p\P\p)P^  +  (Pl  •  P2P\p  -  p\P2p)P^ 

+  [p?P2-(Pi-P2)']^p}  ,  (3.58) 

and 

5"  fe'I^J.p]  =  -8P53P^  .  (3.59) 

Thus,  the  function  ^3  must  satisfy  the  differential  equation 

Ulvz  =  -2e53  .  (3.60) 

3-3   The  Model  Equation  in  the  Asymptotic  Region 

To  obtain  the  solutions  to  the  differential  equations  in  the  asymptotic  region, 

we  must  be  very  careful  how  to  approximate  the  functions  R,Ri,...,S^Si, 

In  the  inner  asymptotic  region,  we  cannot  neglect  those  terms  multipKed  by  the 

electron  propagator  functions  Aj  and  A2  unless  they  are  added  (substracted)  from 

terms  that  are  clearly  larger.  For  example  we  can  approximate  A1A2  +pi  ■p2  ~  Pi  ■p2 

but  cannot  neglect  the  term  (Ai  +  A2)v  against  pi  ■  p2  u.  If  we  do  so,  it  will  not  be 

possible  to  join  the  asymptotic  solutions  to  the  perturbation  solutions.  These  terms 

are  not  negligible  since,  in  the  inner  asymptotic  region,  the  form  of  the  electron 

functions  is 

Ai  ~  m(     xz/m'^)~'^ 

(3.61) 
A2  ~  m{xz~'^ /m^y^ 
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where  by  Eq.(2.38),  r]  =  3e/4  ?a  1.75  x  10"^  Although  the  functions  Ai  and  A2 
approach  zero  for  large  values  of  x,  one  has  go  to  an  extremely  high  momentum 
to  see  an  appreciable  difference.  This  problem  does  not  arise  in  the  first  five  since 
all  terms  containing  the  electron  propagator  functions  are  negligible  compare  to 
those  retained.  In  the  inner  asymptotic  region  the  following  form  for  the  functions 
i?i,  i?2?  S,  S3  are  correct: 

Ri  ~  [v  +  pIu2  -  {pj-pi  -^2)^3]  /D  ,  (3.62) 

R2  ~  [v  +  plui  -  (pI  +  pi  •P2>3]  /D  ,  (3.63) 


r 


s 


I.D  (3.64) 


-Pi  ■  P2V  +  [pIpI  -  {pi  •  P2f] vz 
Sz  r^  [v  +  Pi  ■  P2V:i]  / D  .  (3.65) 

The  remaining  functions  have  the  form 

i?  ~  [pi  •  P2U  +  (Aap2  +  A2P1)  ■  (uiP2  +  U2P1)  +  {Ai  +  A2)v 

+  {PI  -  Pi  ■  P2)v2  -  {pI  -  Pi  ■  P2)vi  -  \pIpI  -  {pi  ■  p2f]uz]  ID  .   (3.66) 

Rz  ~    [-U  -  A1U2  -  A2U1  +pi  ■  p2Uz 

+  {v2-vi)  +  {Ai+A2)vz\lD.  (3.67) 

Si  ~  [-Aiv  -  pIv2  -  {A2PI  +  Aipi  •  P2)vz\  ID  .  (3.68) 

52  ~  [  A2V-plvi+{AiPl  +  A2pi-p2)vz\lD.  (3.69) 

At  this  point,  it  is  important  to  realize  that  in  the  asymptotic  region  the 
system  of  nine  coupled  differential  equations  has  decoupled  into  two  set  of  equations. 
Notice  that  Eqs.(3.62)  to  (3.65)  are  coupled  together  but  not  to  the  remaining  four. 
Therefore,  it  is  possible  to  attempt  to  solve  this  set  of  five  equations  independently 
from  the  rest.  This  is  the  approach  we  will  follow. 

It  is  possible  to  further  simplify  the  differential  equations.  To  this  end  an 
explicit  form  for  the  coefficient  of  the  five  functions  w,  Wi,  ^2,^3,  e  will  be  used. 
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From  Eqs.(3.4),  (3.5)  and  (2.71)  it  follows  that  the  denominator  D  that  appears 
in  the  right-hand  side  of  all  the  differential  equations  can  be  approximated  by  the 
expression 

D^x^  .  (3.70) 

This  approximation  is  valid  in  both  the  inner  and  outer  asymptotic  regions  provided 
that  {x/m?)  ;>  y.  If  we  realize  that  in  the  asymptotic  region,  p^  ~  pi  •  p2  ~  a;y, 
then  the  differential  equations  are  somewhat  simplified.  For  convenience,  all  nine 
asymptotic  forms  of  the  differential  equations  are  presented  in  Table  3-1. 


Table  3-1    Asymptotic  Forms  of  Differential  Equations 
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Olv  -  2h{v)  +  2(ui  +  U2)  =  2e(S  -  e) 
Ulu,  +      2zf_{u,  +  U2)  =  2e[R,  -      z/_(e)] 
Ulu2  +  2z-7+(ui  +  U2)  =  2e[R2  -  z-^U{e)] 
Ulv,  =  -2eSs 
Die  =  5(i?i  +R2)+Pi-nR2+p2-  Di?!  +  h(S) 
Ulu  =  -3eR 


(3.71) 
(3.72) 
(3.73) 
(3.74) 
(3.75) 
(3.76) 


Ulv^  -  ^h{v,)  -      2zf_[p,  •nz;2]  -      2zf_\p2-Uv,] 

-  2zf_{v^+v2)-6u3-2pi-nu3^2€Si  (3.77) 
Olv2  -  4h(v2)  -  2z~'U\Pi  ■  0^2]  -  2z-'U[p2  ■  Dv,] 

-  2z-^f+{vi  +  V2)  +  6M3  +  2p2  ■  Bus  =  2e  ^2  (3.78) 
Ulus  -  10h{u3)  -  2p  -Uhim)  +  h[h{v2  -  v^)]  =-ei23  .      (3.79) 


where 


R  -^  ^[xy  u  +  (A1P2  +  A2P1)  ■  (uip2  +  U2P1)  +  (Ai  +  A2)v 
+  x^/y^-l(vi  +  V2)  +  x\y^  -  l)t/,3] 


Ri 

/-^w' 

i  r 
x^  [ 

i?2 

r\^ 

1    r 

X2    [ 

R^ 


V  +   xz     U2  —  x\/y'^  —  l  U3 

V  +  xz~^  ui  +  xyy"^  —  \vj, 
3  ~  —^  \-u  -  A1U2-  A2U1  +xy  uz  +  {v2  -  vi)  +  (Ai  +  ^2)173] 


5  ~  --  [-y  V  +  a:(y^  -  1)  V3] 
Si  ~  —  [-AlU  -  XZV2  -  x{zA2  +yAi)v3] 
5*2  ~  ^  [    A2V  —  xz~^vi  +  x(2~^  Ai  +  y  ^2)^3] 


[v  +  xy  vz] 


(3.80) 
(3.81) 

(3.82) 

(3.83) 

(3.84) 

(3.85) 

(3.86) 

(3.87) 


CHAPTER  IV 
ASYMPTOTIC  SOLUTIONS  TO  THE  MODEL  EQUATION 

We  seek  solutions  to  the  differential  equations  given  in  Eqs.(3.61)  to  (3.69) 
which  join  smoothly  with  the  perturbation  solutions.  In  order  to  do  this,  we  use 
as  a  guide,  the  asymptotic  form  of  the  perturbation  solutions  found  in  Table  C-2 
of  Appendix  C.  As  can  be  seen  from  this  table,  all  perturbative  solutions  can  be 
written  in  the  form  of  a  power  law  in  the  variable  x.  With  this  in  mind,  we  make 
the  hypothesis  that  the  asymptotic  solutions  to  the  differential  equations  are  of  the 
form 

W^  -  x'^giy)p^  ,  (4.1) 

where,  at  this  stage,  we  will  use  a  bar  over  functions  to  denote  functions  that 
are  only  y-dependent.  The  justification  for  this  assumption  lies  in  the  fact  that 
the  asymptotic  expression  for  the  perturbation  solutions  are  of  this  form  and  they 
satisfy  the  model  equation  in  this  region  (within  .5%)  with  the  exception  of  the 
^1,^2  and  U3  functions.  The  latter  are  good  only  within  3%.  When  we  apply  the 
D'Alambertian  operator  to  such  functions,  it  is  clear  that 

4(«  +  2)r,   2      ^,dg 


[(y2-l)^-ay^]|xV.  (4.2) 


X 

It  is  evident  from  this  equation  that,  if  g{y)  is  a  smooth  function  of  the  variable  y, 
then  for  large  values  of  a;,  it  is  possible  to  neglect  the  first  term  since  it  is  one  power 
of  X  smaller  than  the  second  one.  This  leads  to  the  approximate  form  of  Eq.(4.2) 

D^VF^~.V{^^^.^}  (4.3) 
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where 


ay 


(4.4) 


To  simplify  the  differential  equations  in  the  asymptotic  region,  it  is  necessary  to 
find  the  form  of  the  operators  pi  •  D,  p2  •  0,  ''i,  ^,  etc.,  when  acting  upon  functions 
of  the  form  given  by  Eq.(4.1).  It  is  easy  to  verify  that  when  these  operators  are 
applied  to  functions  of  the  form  $  =  x°'g(y),  the  following  relations  are  valid: 


pi  ■n^  =  x'^ 


zDr,  +  a 


Z     ^£>a  +  « 


y'Da  +  a 


\/y2_ii)„ 


h{^)  =  2x''-^D^g  , 


/($)  =  -2x 


um  =  - 


a-l 


V?'^ 


yDa  -a 


-a-l 


/-(*) = 


V?/'-l 


{Da  -  Oiz)g  , 


{Da  -az     )g  . 


(4.5) 
(4.6) 
(4.7) 
(4.8) 
(4.9) 
(4.10) 

(4.11) 

(4.12) 


Also,  a  combined  application  of  the  h  and  /  operators  upon  these  type  of  functions 
satisfies  the  following  relations: 


h[m]^i[h{^)\ , 

/i[p-n$]  =2A($)+p-n/i($) , 


(4.13) 
(4.14) 
(4.15) 


where 


h\k  .  n$]  =  2x'^-'^y^-l[Dl  -  y  ba\g  , 


(4.16) 


h[p-n^]  =2:r"-i 


yDl+[a-(y'-l)]D, 


i[k  ■  n$]  =  2x^-'  [y  Dl  -  (y2  +  «)^„]  g 


l[H^)]  =  -4- 


X 


a-2 


y  Di-\{y'-l)  +  a  D 
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(4.17) 
(4.18) 
(4.19) 


and  the  operator  £)^  is  defined  by  the  equation 

Dl  -  iv'  -  1)'!^  -  2(a  -  l)y(y2  -l)±  +  a[l  +  {a-  l)y']    .  (4.20) 

Using  Table  C-2  of  Appendix  C,  it  is  reasonable  to  assume  the  following  forms  for 
the  functions  u,  u^,  U2,  v^  and  e: 


Ui   =  —U2  —  X     U\ 


and 


V^^X     Ug    , 


e  =  a;"e 


(4.21) 
(4.22) 
(4.23) 

(4.24) 


where  c^  is  an  arbitrary  constant.  Furthermore,  the  asymptotic  behavior  of  the 
u-function  (see  Appendix  D)  leads  to  the  conclusion  that  the  value  of  a  must  be 
equal  to  -1.  Other  solutions  are  possible  with  a  ^  -I  but  they  proved  to  vanish 
at  large  values  of  the  variable  x.  Such  solutions  are  not  favored  since  they  do 
not  have  an  asymptotic  behavior  consistent  with  Appendix  D  and  they  do  not  join 
smoothly  with  the  perturbation  formulas  found  in  Eqs.(C.41)  to  (C.46).  To  find  the 
asymptotic  solutions  to  equations  (3.71)  to  (3.75)  that  join  the  perturbation  forms, 
we  set  all  functions  in  the  right-hand  side  of  the  differential  equations  equal  to  zero 
except  for  the  v  and  e-functions.  Substituting  Eqs.(4.21)  to  (4.24)  into  Eqs.(3.71) 
to  (3.75)  and  setting  a  —  —1  gives 


4[L>_i  +  y]v  =  -2e[yv  +  e]  , 


(4.25) 
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4D_i«i  =  2e[v  +  z{e  +  ^/y^  -  1  e')]  ,  (4.26) 

4I)-iU2  =  H^  +  ^~^(e  -  Vy'  -  1  e')]  ,  (4.27) 

4I)_ii;3  =  -2ev  ,  (4.28) 

and 

4I)_ie  =  10?;+|     z(I)_i  -  y)  -  2}  v^\^z-\b_-,  -  y)  -  2}  y+2D_i(-yu)  (4.29) 
where 

^-i=(i-y')^-y.  (4.30) 

The  primes  denote  derivatives  with  respect  to  the  variable  y. 

Notice  that  with  the  choice,  a  =  —1,  all  five  functions  are  in  agreement 
with  the  X- dependency  of  the  perturbation  solutions  and  in  effect  have  reduced 
the  problem  to  that  of  finding  the  solutions  to  Eqs.(4.25)  to  (4.30)  which  are  only 
y-dependent.  To  solve  these  equations,  we  start  by  adding  Eqs.(4.26)  and  (4.27). 
This  gives  the  result 

2r>_i(fxi+«2)  =  e(^--C>_ie)  .  (4.31) 

Since  we  assume  u\  =  —U2  in  obtaining  these  differential  equations,  we  conclude 
that,  for  consistency, 

D-ie  =  v.  (4.32) 

Applying  D_i  to  Eq.(4.25)  and  substituting  the  previous  expression  gives 

4Z>_i  [£)_i  +  y]v  =  -2e[b{yv)  +  v]  .  (4.33) 

Using  the  identity 

by  =  yb-{y''  -I)  ,  (4.34) 
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and  moving  all  the  terms  to  the  left  side  of  in  Eq.(4.33)  yields  the  result, 

4£'_i [D_i  +  y]v  +  2e[yD_iv  -  {y^  -  l)v]  +  2eu  =  0  .  (4.35) 

In  terms  of  the  operator  D^^  obtained  from  Eq.(4.20)  by  letting  a  =  -1, 

Dl^v  +  (1  +  le)yD_,v  +  [(1  -  le)(y^  -  1)  +  ^e]v  =  0  .  (4.36) 

In  terms  of  the  derivatives  of  v  with  respect  to  y,  Eq.(4.36)  can  be  written  as 

iy^-l)v"  +  {3~le)yv' -ev  =  0  .  (4.37) 

It  is  possible  to  transform  Eq.(4.37)  with  the  substitution, 

v  =  {y^  -  lyv  ,  (4.38) 

v'  ^  {y^  -  lyv'  +  2Ty{y^  -  ly-'  V  ,  (4.39) 

and 

z;"  =  (y2-l)-^"  +  4ry(y2_i)-i^'  +  {2r(y2_i)-i+4T(r-l)j/2(y2_i)r-2|=_ 

(4.40) 
Substituting  Eqs.(4.38)  to  (4.40)  into  Eq.(4.37)  gives 

(y^  -  l)^|(y2  -  1)^"  +  [At  +  (3  -  '^e)]yv' 

+  [2t  +  [4r(r  -  1)  +  2r(3  -  ^e)  ]y\y^  -  ly'  -  e]  ^|  =  0  .       (4.41) 

To  simplify  this  equation,  we  select  the  value  of  the  constant  r  to  satisfy  the  con- 
dition 

4r  +  (3-ie)  =  2  (4.42) 

or  solving  for  r 

r  =  i(-l  +  |6)  .  (4.43) 
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Making  these  substitutions  in  Eq.(4.41)  gives 

(y^-l)v"+2yv'  +  {2T  +  [4t(t  -  1)  -  2t(At  -  2)]y\y^  -  1)"^  _  e}^  =  0  ,  (4.44) 
or  after  some  simplifications 

(2/2  -l)v"  +  2yv'  -  {(4r2  -  2r  +  e)  +      2^_  i)}  ^  =  ^  •  (4-45) 

If  we  define  the  constant  /j.  such  that  //  —  ±2r  and  the  constant  u  by  the  relation 

u{u  +  1)  =  4r2  -  2r  +  e  ,  (4.46) 

then  Eq.(4.45)  can  be  written  in  the  form 

(y^  -l)v"  +  2yv'  -  {p(u  +  1)  +  .  /    ^    1  y  -  0  ,  (4.47) 

(y    -  1)  ^        ^ 

where 

/i  =  ±|(-1  +  |e)         and         z^  =  |(1  +  |e)  .  (4.48) 

The  general  solution  to  Eq.(4.47)  can  be  expressed  in  terms  of  the  Associated  Leg- 
endre  functions  P^^(z)(See  Appendix  E).^^'^^  These,  in  turn,  can  be  represented  by 
the  hypergeometric  series  F(a,  b,  c;  ij)  by  the  identity, 

^"^'^  =  ''^'r(7-,")'''^^^^'' - 1^, I  - 1- - 1/^,  1  - /^;  1  - y~') ,    (4.49) 

vaHd  for  y  >  1.  For  computational  purposes,  however,  this  expresion  is  not  conve- 
nient due  to  the  slow  convergence  of  the  hypergeometric  series  when  its  argument 
is  near  1.  For  large  values  of  y  the  following  expression  gives  a  rapid  convergence 
of  the  series: 

my)  =      r    r    ' — W~'-\v'  - 1)-''' 
2"    r(i  +  zy)      ^     ^ 

a/tt  i  (1  +  1/  -  ^) 

X  F{-\iX  -\y^\-\y-  i^,  i  _  y.  y-^)   .  (4.50) 
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The  u-function  can  then  be  written  in  the  form 

v  =  c,{y^~iy^^Pi:{y).  (4.51) 

where  /i  and  v  are  defined  in  Eq.(4.48).  We  now  turn  to  the  problem  of  solving 
Eq.(4.32)  with  the  solution  to  the  u-function  differential  equation  in  the  right-hand 
side.  From  Eq.(4.32)  and  the  expression  for  D-i  given  in  Eq.(4.30),  the  differential 
equation  that  the  e-function  must  satisfy  becomes 

(r/  -  1)  e '  +  y  e  =  -c,{y^  -  l)'''P,''(y).  (4.52) 

Let  assume  that  this  function  has  the  same  form  as  the  u-function,  i.e., 

e  =  -ci(y2-l)^P;(y)  (4.53) 

and 

e'  =  -2Ay(y2  -  \f-'P^{y)  -  (y^  -  l)P^'{y)  .  (4.54) 

Substituting  these  two  equations  into  the  preceding  one  gives 

(y^  -  1)^|(2A  +  l)yP;(y)  +  (y^  -  l)P;'(y)}  =  (y^  -  l)'''P,^(y)  ,         (4.55) 
from  which  it  follows  that  A  =  |/^.  With  this  substitution,  this  equation  becomes 

{^i  +  l)y  P^{y)  +  (y2  -  1)P«  '(y)  =  Pl^{y)  .  (4.56) 

Using  the  identity, 

(y^  -  l)PI  '(y)  ^  /?  y  P^{y)  -  (a  +  p)P^_,{y)  ,  (4.57) 

yields  the  equation 

(/.  +  /?  +  l)y  P^{y)  -{a  +  /3)P^_,(y)  =  Pj^iy)  .  (4.58) 
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If  we  set  a  =  ;U,  and  /?  =  -(1  +  /i)  then  it  follows  that  the  first  term  in  this  equation 
vanishes  and  the  conditions, 

-(«  +  /?)  =  !         and     ^-l  =  -(z.  +  i),  (4.59) 

are  satisfied.  It  is  evident  from  Eq.(E.17)  that  the  e-function  is  given  by  the  equation 

e  =  -C^y"  -■l)^''P!\i+,,)(y)  =  -ci(y2  -  l)'^P;(y)  .  (4.60) 

We  now  turn  to  the  problem  of  solving  Eq.(4.26)  to  find  the  function  u^.  To  do 
this,  we  will  transform  this  equation  with  the  help  of  Eq.(4.12).  From  this  equation, 
(with  a  —  —1),  it  follows  that 

^(e  +  v'y'-le')  = rir (D_i+.g-^)e.  (4.61) 

V  y^  -  1 

Substituting  Eq.(4.32)  into  this  equation  and  simplifying  gives 

z(e  + Vy'-le')  =  — 7=^-  (4.62) 

With  this  result,  Eq.(4.26)  yields  the  equation 

"      _                  2e 
4D_i«i  = {xjv  +  e]  ,  (4.63) 

vy  -1 

or  in  terms  of  the  first  derivative  of  ui, 

{y" -\)u[-Vyu^^-e—=^\yv^-t\.  (4.64) 

^    \Jy^  -  1 

To  solve  this  equation  we  make  the  substitution, 

^1  =  ci(2/2  -  1)^(^-1)^1  (4.65) 

and 

u[  =  ci(/i  -  l)y^i(y2  -  l)i(''-^)  +  (y2  -  l)^('^-i)(y2  -  1)^^  '  .  (4.66) 
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Then  the  left-hand  side  of  the  Eq.  (4.64)  becomes 

{y^  -  l)u{  +  yui  =  (y2  _  i)i(/'-i)  |^  y=^  +  (^^2  _  ^^  -^  '1    _  ^^_g^^ 

The  right-hand  side  of  Eq.(4.62)  can  be  written  in  the  form, 

where  we  have  substituted  the  expressions  for  v  and  e  found  in  Eqs.(4.51)  and 
(4.60).  Equating  the  expressions  given  in  Eqs.(4.67)  and  (4.68)  yields  the  equation 

(y^  -1)Z^'  +  f.t  yu,  =  ^  [yP;:(y)  -  P;(y)]  (4.69) 

where  ^  and  u  are  defined  in  Eq.(4.48).  Notice  also,  that  from  Eq.(4.48),  u-l  =  fi 
and  fj.  +  u^  ^e.  Therefore  Eq.(4.69)  can  be  rewritten  in  the  form, 

(y'  -l)u,'  +  fiyZ,=(^i  +  v)  [yPJiiy)  -  P^^{y)\   .  (4.70) 

It  is  possible  to  transform  the  right-hand  side  of  the  last  equation  by  using  several 
relations  satisfied  by  the  Legendre  functions  (see  Appendix  E;  Eqs.(i;.43)  to  (E.51)). 
It  is  straightforward  to  show  that 

(y2  -  1)^1 '  +  ^  j/^^  =  (^  _  ^  +  l)P.\i(y)  -  2yP,''(y) 

=  {y'-\)P^{y)'  -^liyP^iy).  (4.71) 

It  is  easy  to  verify  that  the  general  solution  to  the  homogeneous  equation  is  given 
by  the  expression 

(^i)/.=C2P;_i(y)  (4.72) 

where  C2  is  an  arbitrary  constant.  It  is  evident  from  Eq.(4.71)  that  a  particular 
solution  to  this  equation  is 

^^)v  =  P^{v)  ■  (4.73) 
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If  we  set  the  constant  C2  =  —1,  then  the  general  solution  to  Eq.(4.71)  is 

^1  =  P.'iy)  -  Pll-M  =  (2a^  +  i)Vy'-ip;-^(y) .  (4.74) 

Using  Eqs.(4.65)  and  (4.66)  we  finally  arrived  at  the  desired  solution 

ui  =  -U2  =  ci(2^  +  l)(y2  _  l)l''p^^-i(j/)  .  (4.75) 

The  i;3 -function  can  be  easily  obtained  if  we  reaHze  that  the  differential  equation 
that  it  must  satisfy  is  very  similar  to  the  differential  equation  satisfied  by  the  e- 
function.  In  fact,  we  notice  that  the  difference  between  Eq.(4.28)  and  Eq.(4.32)  is 
a  constant  factor.  With  this  in  mind,  the  function  U3  can  be  expressed  in  terms  of 
the  e-function  in  the  form, 

v^^{li  +  v)e  =  -{^x  +  z/)ci(y2  _  l)l'^P;(y)  .  (4.76) 

The  arbitrary  constant  ci  needed  to  completely  define  the  solutions  can  be 
determined  by  using  the  requirement  that  the  u-function  must  approach  the  value 

^  -  (1  -  ^e)  =  (1  -  \ri)  (4.77) 

as  a;  — >■  00  and  y  ^-  1.  This  result  was  derived  in  Appendix  D.  From  Eq.(4.51)  and 
(E.57)  we  have  that 

lim  Pl^iy)  =        ^'       (y^  -  1)-^'^  (4.78) 


and  therefore 


2^ 
lim^'  =  Ci=;T- r.  (4.79) 


Equating  the  previous  expression  with  Eq.(4.75)  yields  the  result 

,,     1  ,r(i-^) 
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We  now  turn  to  the  problem  of  solving  Eq.(3.76)  and  (3.80)  for  the  function  u. 
Its  is  possible  to  simplify  the  function  R  given  in  Eq.(3.80)  in  the  following  manner. 
If  we  expand  the  scalar  product  involving  the  functions  Uj  and  U2  and  substitute 
the  condition  that  U2  =  -wi,  this  equation  becomes 

R=  ^Ixyu  +  ui  [Aipj  +  A2P1  ■  P2  -  ^iPi  •  P2  -  A2PI] 

+  (Ai  +A2)v  +  x^/y^  -  1  (ui  +  V2)  +  x^(y^  -  l)u3  1  .  (4.81) 

Using  the  asymptotic  forms  oipj,  pj,  and  pi  ■  p2  given  in  Eqs.(3.4)  and  (3.5)  yields 
the  result 

^  =  ^1^^  "  +  ("  -  Vy^-'^Ui){Ai  +  A2)v  +  x^y^-l{vi  +  V2)  +  x\y^  -  1)^3 1 

(4.82) 
The  terms  involving  the  functions  v  and  Ui  can  be  simplified  further.  Substituting 
the  expressions  for  these  functions  found  in  Eqs.(4.51)  and  (4.75)  gives 

V  -  Vy^-'^u,  =  c,{y^  -  1)^^  ]^PI^{y)  -  (2fi  +  l)v/y2-lP;-l(y)]    .         (4.83) 

Using  equation  (E.46)  found  in  Appendix  E,  it  follows  that  the  previous  equation 
can  be  written  as 

V  -  V^^^u,  =  c,(y'  -  1)^/^  [Pl^iy)  +  p;:_^(y)  _  p;^^(y)]    .  (4.84) 

If  we  notice  that  fj,  +  l  =  u,  this  expression  becomes 

V  -  Vv^-luj  =  ci(y2  _  l)i''P^^(y)  .  (4.85) 

Furthermore,  using  Eq.(£'.41)  it  follows  that 

v-^y-'-lu,  =ci2Vr(l-/i)  .  (4.86) 


51 

or,  using  Eq.(4.80) 

y  -  Vy^  -  1  Ml  =  (1  -  -,7)  .  (4.87) 

If  we  substitute  this  expression  in  Eq.(4.83)  and  set  all  other  functions  equal  to 
zero,  we  arrive  at  the  equation 

Ulu  =  -3e(l  -  ^77)(Ai  +  A2)/x2  .  (4.88) 

Using  the  same  approach  used  in  solving  the  first  five  differential  equtions,  we 
assume  that  the  u-function  can  be  represented  by  the  form 

u  =  x'^u{y)  .  (4.89) 

Substituting  the  previous  form  into  equation  (4.88)  and  using  the  asymptotic  form 
of  the  D'alambertian  operator  as  well  as  the  asymptotic  forms  of  the  functions  Ai 
and  A2  gives 


x^-H{/3  +  2)b^u  =  -3e(l  -  \r])  x-'^-'>{z^  +  ^-'')  . 


2 


(4.90) 


Matching  the  x-dependency  on  both  sides  of  Eq.(4.90)  implies  that  the  parameter 
^  must  be  defined  by  the  formula 

/9  =  -(l+77).  (4.91) 

Substituting  the  form  of  the  I)  given  in  Eq.(4.4)  yields  the  equation 

{x/  -  l)ii '  +  (1  +  r])yu  =  A(z''  +  ^-'')  ,  (4.92) 

where 
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To  solve  this  equation,  we  notice  that,  for  values  of  y  >  10,  the  approximation  z  ^  2y 
is  very  good.  This  leads  to  a  relative  error  of  less  than  .2%.  Since  the  perturbation 
solutions  are  good  in  this  region,  we  will  attempt  to  solve  the  approximate  equation 


y^u'  +  {l+rj)yu  =  X[{2yy  +  (2y)-^]    . 
In  order  to  solve  this  equation,  we  make  the  change  of  variables 
^  =  {2y)\      d^^2rj^^'-^\      and     y  =  ^^^  ■ 


(4.94) 


(4.95) 


Subtituting  these  expressions  in  Eq.(4.94)  and  simpHfying  yields  the  differential 
equation 

u*  +  {l  +  Ih)  J  =  2{X/ri)rHl  +  r']  ,  (4.96) 

where  the  symbol  u*  stands  for  the  derivative  of  u  with  respect  to  ^.  This  equation 
can  be  solved  by  introducing  the  integrating  factor 


1(0 -exp|(i +  1/^7)  III  =e(^+i) 


(4.97) 


The  general  solution  can  be  written  in  terms  of  the  integrating  factor  in  the  form 

T  1         r^ 

^  =  7(^  +  (2V^)j(^y  {t-Hi+t-']}mdt,         (4.98) 

where  d  is  an  arbitrary  constant.    The  integration  in  the  previous  equation  is  ele- 
mentary an  yields  the  result 


u  =  dr^'  +  ^^  +  {2X/7j)r^^+^^ 


ini+-e 


(4.99) 


In  terms  of  the  variable  y,  the  solution  can  be  written  as 


^  =  ^  H^y)'"  +  {^/ri){2yr  +  2A(2y)-''  ln(2y)}    . 


(4.100) 


To  determine  the  constant  d,  we  require  that  this  solution  join  smoothly  with  the 
asymptotic  form  of  the  perturbation  solution  given  in  Eq.(C46)  of  Appendix  C. 
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In  order  to  join  the  solution  just  found,  we  notice  that  for  values  y  not  too  large 
the  following  expansions  are  valid  due  to  the  small  magnitude  of  the  constant  rj 
{■q  ^  1.75  X  10-3): 

{2yy     =exp{    ?7ln(2y)}  =  l  +  ?7ln(2y)  +  ---   , 

(4.101) 
(2y)-''  =  exp{~r}ln{2y)}  =  1  -  rj\n{2y)  +  ■■■   . 

If  we  select  the  value  of  the  constant  d  to  be  equal  to  -X/rj  ,  then  Eq.(4.97)  becomes 
^  =  ^{^[(2y)''-(2j/)-'']  +  2(2y)-''ln(2y)|   ,  (4.102) 

or  because  of  the  expansion  given  in  Eq.(4.101), 

2A 


w  ->  —  ln(2y)[l  -  T]  ln(2y)]  .  (4.103) 

Substituting  the  expression  for  A  given  in  Eq.(4.90)  gives 

3   ln(2y)  ^^,  „, 

n  -  l^e-Y^    +   Oie')  ,  (4.104) 

^      y 

which  has  the  same  form  as  that  given  in  Eq.(C46).  Therefore,  the  solution  is 
found  to  be  given  by  the  expression 

^  "  ^^^  "  « '^  ^^^^^^  ~  ^^y^~'  +  2r?(2y)-''  ln(2y)}   .  (4.105) 

It  is  not  difficult  to  construct  a  "solution"  that  is  accurate  near  y  =  1.  From  the 
form  of  the  perturbation  solution  near  y  =  1  and  the  fact  that  we  can  replace  all 
terms  involving  z  by  (2y),  it  is  reasonable  to  assume  that  this  replacement  wiU 
improve  the  solution  in  this  region.  With  this  argument  the  improved  solution  is 
given  by  the  expression 

^  =  ^(1  -  1^)  K^"  -  ^~')  +  2^^""  Inz}  /\/^^^  .  (4.106) 
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For  reference  purpose,  the  asymptotic  solutions  are  presented  in  the  Table  4-1. 
Numerical  values  for  the  functions  u,  e,  ui,  v^  and  u  are  presented  in  Table  4- 
2.  The  variable  x  inside  the  overlapping  region  was  chosen  to  be  10^°.  Their 
derivatives  are  also  presented.  The  values  of  the  perturbations  solutions  are  also 
shown  for  comparison. 


Table  4-1    Asymptotic  Solutions  to  the  Model  Equation 
Notation: 

//=|(-1  +  Je);  ^=i(l  +  ie); 


,^-e 


v^{i-  lv)iy'  -  it^'P.'iy) .  (4.107) 

ui  =  -U2  =  (1  -  ^77)(2^  +  l)(y2  _  l)l/^p;-i(y)/a;  . 

^3  =  -ie(i  -  lv){y'  -  l)^''P;(y)/x  .  (4.110) 

u  =  ^(1  -  le)x-^'+^^  {(^^  -  z-^)  +  2^z-^  \nz}  /Vy^^  ■  (4.111) 


(4.108) 
(4.109) 
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Table  4-2    Perturbation  and  Asymptotic  Solutions  in  the  Overlapping  Region 


1**  derivative 
2"''  derivative 


Perturbation  Solutions  (x  —  lO'^*'  and  y  =  1.1) 


9.957E-01 

7.301E-04 

-4.155E-04 


-9.679E-01 

3.083E-01 

-2.352E-01 


1.640E-04 

7.341E-04 

-4.877E-03 


^^3 

1.124E-03 

-3.579E-04 

2.731E-04 


3.372E-03 

-1.074E-03 

8.193E-Q4 


1*'  derivative 
2"^*  derivative 


Asymptotic  Solutions  (x  =  10^°  and  y  —  1.1) 


9.992E-01 

7.298E-04 

-4.149E-04 


-9.671E-01 

3.078E-01 

-2.34SE-01 


1.640E-04 

7.338E-04 

-4.S75E-03 


^3 

1.123E-03 

-3.574E-04 

2.726E-04 


3.373E-03 
-1.081E-03 

8.574E-04 


Perturbation  Solutions  {x  —  10^"  and  y  =  10) 


P*  derivative 
2"^^  derivative 

V                        e 
9.979E-01       -3.008E-01 
1.138E-04        2.029E-02 
-l.lOlE-05       -3.109E-03 

Asymptotic  Solutions  (x  = 

2.344E-04 

-1.224E-05 

1.390E-06 

10^°  and  y  = 

V3 

3.493E-04 

-2.355E-05 

3.609E-06 

10) 

u 

1.048E-Q3 

-7.066E-05 

1.083E-05 

1**  derivative 
2'^^  derivative 

V 

l.OOlE+00 

1.139E-04 

-1.102E-05 

e 

-3.008E-01 

2.027E-02 

-3.105E-03 

2.345E-04 

-1.223E-05 

1.388E-06 

^^3 

3.493E-04 

-2.354E-05 

3.606E-06 

u 

1.046E-03 

-7.063E-05 

1.083E-05 

Perturbation  Solutions  (x  —  10      and  y  =  100) 


P*  derivative 
2"    derivative 

V 

l.OOlE+00 

1.161E-05 

-1.160E-07 

Asymptotic 

e 

-5.299E-02 

4.299E-04 

-7.599E-06 

Solutions  (x  = 

Si 

4.991E-05 

-3.831E-07 

6.504E-09 

lO^o  and  y  = 

V3 

6.152E-05 

-4.992E-07 

8.823E-09 

100) 

u 

1.846E-04 

-1.497E-06 

2.647E-08 

I''*  derivative 
2^^  derivative 

V 

1.004E+00 

1.165E-05 

-1.163E-07 

e 

-5.305E-02 

4.301E-04 

-7.601E-06 

5.000E-05 

-3.S35E-07 

6.508E-09 

6.160E-05 

-4.994E-07 

8.825E-09 

1.839E-04 

-1.493E-06 

2.641E-08 
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Table  4-1    (Continued) 

Perturbation  Solutions  (x  =  lO'^^  and  y  =  1000) 

V  e                       ui                        vs  u 
1.003E+00       -7.601E-03         7.664E-06        8.825E-06  2.648E-05 

1^*  derivative           1.161E-06         6.601E-06       -6.503E-09       -7.664E-09  -2.299E-08 

2"^^  derivative        -1.161E-09       -1.220E-08         1.185E-11         1.417E-11  4.251E-11 

Asymptotic  Solutions  (x  =  10^"  and  y  =  1000) 

V  e                       ui                        V3  u 
1.007E+00       -7.620E-03         7.687E-06        8.848E-06  2.633E-05 

1^*  derivative          1.169E-06        6.613E-06       -6.518E-09       -7.679E-09  -2.288E-08 

2"^^  derivative        -1.168E-09       -1.222E-08         1.187E-11         1.419E-11  4.232E-11 


CHAPTER  V 
THE  CONCLUSION 

A  solution  for  the  vertex  amplitude  has  been  found  to  the  Schwinger-Dyson 
equations  based  on  an  approximation  scheme  which  is  characterized  by  the  follow- 
ing: 

(1)  the  photon  propagator  is  approximated  by  its  form  near  the  mass  shell, 

(2)  the  infinite  hierarchy  of  the  vertex  is  cut  off  at  the  second  order  in  the 
couphng  constant  and  the  remainder  is  approximated  by  Green's  generahzation  of 
the  Ward  Identity  for  higher  order  contributions. 

The  gauge  is  chosen  to  be  that  found  necessary  to  obtain  a  finite  solution  to 
the  electron  propagator  equation^"  with  a  vanishing  bare  mass  for  the  electron.  In 
the  above  reference,  the  electron  propagator  is  found  to  be  given  by 

S-\p)^-A{p')+B{p')P  (5.1) 

where 


A{p')  = 


-2  3a(m^-p^)/47rp^ 


^-1 


(5.2) 

Bip")  =  1  .  (5.3) 

In  the  asymptotic  region,  a  simpler  form  of  the  A-function  was  used,  namely 

A(p2)  ~  {fjm^)-'^     r]  ^  1.75  x  10"^  .  (5.4) 

A  simplified  approximation  of  the  vertex  equation  (the  model  equation)  was 
found.     This  tensor  and  matrix  equation  can  be  reduced  to  a  set  of  8  coupled 
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differential  equations  of  third  order  in  two  variables.  The  order  of  the  equations 
was  reduced  by  the  introduction  of  an  additional  function  e.  The  set  of  8  basic  scalar 
functions  is  needed  to  define  the  vertex  completely  as  required  by  the  transformation 
properties  of  the  vertex  amplitude.  In  general,  the  transverse  part  of  the  vertex 
amplitude  can  be  written  in  the  form 

r  iPuP2)=J_^v{pl,pl,pj  ■P2) 
+  P^u{pI,pI,Pi  ■P2) 

+  Ap^Mp1^pI,pi  •P2) 

+  ^2P^U2(pI,pI,Pi  ■P2) 

+  ie^'""''l^PiuP2pV^{pl,pl,p^  ■P2)  .  (5.5) 

The  longitudinal  part  of  the  vertex  is  given  by  Ward's  identity  (see  Eq.(2.27)). 

The  analytical  approximate  solution  to  this  simpler  equation  has  been  found 
in  the  vicinity  of  the  four  momenta  square  equal  to  the  square  of  the  experimental 
mass.^^'^^  The  validity  of  these  solutions  was  tested  by  direct  substitution  into 
the  differential  equations  (see  Eqs.(3.71)  to  (3.79)  with  Eqs.(2.85)  to  (2.92)).  A 
measurement  of  their  accuracy  is  given  by  the  difference  between  the  left  and  right- 
hand  side  of  the  differential  equations  divided  by  their  average.  These  solutions 
were  found  to  be  vaHd  within  .5%  for  6  of  the  equations  (u,  Ui,  U2,  ^3,  e,  and  u) 
and  a  much  higher  error  (3%)  for  the  remaining  functions  (ui,  ^2,  1*3).  The  region 
of  vaHdity  of  these  solutions  was  found  to  be  a;  >  1  and  1  <  y  <  10^. 
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The  introduction  of  the  variables, 


^! 


X  =  ^pIpI         and         z  =  ^Jpl/pl  ,  (5.6) 

made  it  possible  to  separate  the  system  of  9  coupled  differential  equations  into  three 
separate  systems  of  five,  three  and  one  differential  equations  in  the  asymptotic 
region  of  large  momenta.  Using  these  variables,  the  solutions  factor  into  functions 
of  X  and  functions  of  z  in  this  region,  and  the  x-dependency  of  all  the  solutions 
can  be  written  in  a  simple  power  law  form.  This  reduces  the  problem  to  finding 
functions  of  only  one  variable. 

The  goal  of  this  work  was  to  find  the  solutions  to  the  set  of  differential  equa- 
tions containing  the  largest  component  of  the  vertex  amplitude,  namely  the  v- 
function.  The  differential  equations  for  the  functions  u,  wi,  U2,  W3,  e,  are  linked 
together  but  do  not  involve  the  other  dependent  variables  (see  Eqs.(3.71)  to  (3.75)  ). 
Solutions  were  found  for  these  five  functions  that  join  smoothly  with  the  asymptotic 
forms  of  the  perturbation  solutions.  The  percent  difference  between  the  perturba- 
tion and  asymptotic  solutions  is  less  than  .1%  in  the  region  x/k'^  >  10-^°  and 
1-1  ^  y  ^  10^  with  a  similar  result  for  their  first  and  second  derivatives.  This 
confirms  the  results'^^  already  found  that  the  perturbation  solutions  gives  a  good 
representation  of  the  vertex  amplitude  in  a  very  large  region  of  the  electron  mo- 
mentum. This  conclusion  follows  from  the  fact  that  the  largest  part  of  the  vertex 
amplitude  (z.e.  v^  )  probably  has  a  maximum  at  the  mass  shell.  The  smalkiess  of 
the  fine  structure  constant  and  the  slowness  of  the  decay  of  the  electron  propagator 
function  A  in  the  asymptotic  region  are  also  responsible  for  this  result. 

These  results  were  encouraging  so  that  an  attempt  was  made  in  solving  the 
ti-function.  The  techniques  and  general  insight  gained  in  the  solution  of  the  v- 
function  proved  useful.  This  attempt  was  successful  and  a  solution  was  found  that 
also  joined  smoothly  to  the  perturbation. 
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The  form  of  these  functions,  except  the  u-function,  vanishes  at  infinity  as 
rapidly  as  0{l/x)  and  therefore  is  neghgible  compared  to  the  u-function.  Also, 
these  functions  are  one  order  in  the  couphng  parameter  smaller  than  the  u-function. 
We  are  lead  to  conclude  that  the  most  important  function  in  defining  the  vertex 
amphtude  is  v.  Although  the  perturbation  solution  gives  a  good  representation  of 
this  function  in  the  intermediate  region,  this  form  may  not  be  useful  in  determining 
the  convergence  of  several  interesting  quantities.  The  analytical  properties  of  the 
vacuum  polarization  integral  and  charge  renormahzation  constant  Z3  are  strongly 
dependent  of  the  analytical  properties  of  the  vertex  amplitude.  The  convergence  of 
this  integral  may  well  depend  on  the  way  the  vertex  amplitude  behaves  at  infinity. 
The  perturbation  expression  given  in  equation  (C.41)  shows  a  logarithmically  diver- 
gent form  for  this  function  as  the  variable  z  growth  without  bounds.  The  present 
work  has  demonstrated  that  the  behavior  is  much  less  rapid  {y  ~  z^^). 

This  result  alone  encourages  the  calculation  of  the  vertex  amphtude  to  higher 
accuracy  to  determine  if  it  is  constant  everywhere  or  even  decreases  for  large  values 
of  the  variables.  To  address  this  question,  a  recalculation  of  the  photon  propagator 
is  needed.  The  form  of  the  solutions  found  in  this  work  may  lead  to  the  calculation 
of  a  correction  to  the  photon  propagator.  This  in  turn  will  provide  a  way  to  calculate 
the  renormahzation  constant  Z^  and  the  bare  charge  of  the  electron. 

It  cannot  be  said  that  the  present  solutions,  described  in  this  work,  consti- 
tute a  complete  resolution  of  the  problems  in  the  theory  of  QED.  However,  this 
method  has  been  successful  in  finding  these  solutions  without  the  unreasonable  in- 
finite quantities  that  plague  renormahzation  theory.  Its  felt  that  within  the  level  of 
approximations  made,  the  present  solutions  give  evidence  that  the  infinite  quantities 
that  occur  in  the  usual  perturbation  calculations  of  the  self-energy  of  the  electron 
and  the  vertex  amplitude  are  not  essential  to  the  theory.  Furthermore,  the  present 
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work  sets  the  foundations  for  the  subsequent  calculation  of  the  photon  propagator 
and  the  closing  of  the  iteration  procedure  of  the  project. 


APPENDIX  A 
DIRAC  MATRICES:  DEFINITION  AND  IDENTITIES 

The  relativistically  invariant  equations  which  in  the  case  of  the  electron  play 
the  same  role  as  the  Maxwell's  equations  do  in  the  case  of  the  photon,  were  obtain 
by  P.A.M.  Dirac^*^  in  1927.  In  order  to  describe  these  equations  in  a  relativisti- 
cally invariant  form,  Dirac  introduced  a  set  of  four-dimensional  matrices  (the  Dirac 
matrices)  given  by 

I°=(J    1]-'     l'={-l    f)    -=1.2,3  (Al) 

where  a^  represents  the  Pauli  spin  matrices.  The  Pauli  spin  matrices  are  two- 
dimensional  matrices  defined  by 

The  Dirac  gamma  matrices  satisfy  the  anticommutator  relation, 

The  multiple  products  of  gamma  matrices  form  a  group  of  16  linearly  independent 
4x4  matrices.  One  possible  representation  of  this  group  is  given  by  the  linear 
combination  of  products  defined  by  the  matrices /,  7   ,7     ,7         where 


and 
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If 

and  I  is  the  4-dimensional  unit  matrix.  We  have  used  the  symbols  [,]  and  {,}  to 
represent  the  commutator  and  anticommutator  respectively.  Another  useful  matrix, 
used  in  defining  the  model  equation  is  defined  as 

7'  =  -ZT^'l'l"  •  (A.6) 

Equation  (A.6)  may  be  used  to  write  Eq.(A.5)  in  the  form 

%up  =  ^■^/^-pc.lsT"  ■  (A.7) 

The  symbol  e^^p^  stands  for  the  4-rank  antisymmetric  tensor  (Levi  Civita)  defined 
by  the  rules 

Sfippa  —  0         unless  ^u,  u,  p,  a  are  all  different  ,  (^-8) 

£1234  =  1  ,  (A.9) 

fivpcr  _  f  +1      for  even  permutation  of  the  indices  1,2,3,4; 

\  — 1      for   odd  permutation  of  the  indices  1,2,3,4;  ^    '     ^ 


and 


^  ixvpcr 


—  1      for  even  permutation  of  the  indices  1,2,3,4; 
+  1      for    odd  permutation  of  the  indices  1,2,3,4. 


(All) 
It  is  simple  to  verify  that  the  gamma  matrices  satisfy  also  the  following  identities: 

l\x  =  4  '  (^-12) 

t\^1x  =  -27„  ,  (A.13) 


-Lfi-LX  J-n 


l\.lx  =  0  '  (^-14) 


iX.plx  =  27,,,  ■  (A.15) 
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By  definition,  the  operator  ^  =  7^ 9^  and  therefore,  using  Eqs.(A.12)  to  (A. 15)  it 
is  possible  to  show  that 

%,  =  -!/ +  25,,  (A16) 

^1..=!,.^  + 2(5.1. -^''7,),  (A17) 

^I..P  =  -1,J  +  2(5.7,,  +  5.7,^  +  5,7,  J  .  (A.18) 


APPENDIX  B 
DERIVATION  OF  GREEN'S  PERTURBATION  SOLUTIONS 

To  obtain  the  perturbation  solutions  to  the  model  equation,  one  sets 
Cj  ~  ^i.  ^C^  ~  C^^  ~  C^^p  ~  0  and  Ai  ~  7I2  ~  m  as  the  zero*''  order  ap- 
proximation in  equations  (2.52)  to  (2.55).  Therefore,  the  D  tensors  become 

D^  -  2mp^/D  ,  (B.l) 

Di=     [{m'-p^-p2y6^^+ip,,+p2,)p^]/D  ,  (B.2) 

D^,  =   m  [(m^  -  pi  .  ;?2)^;  -  (pi^  -  P2^)^i]  /D  ,  (B.3) 

D^up  =     ^  [{PliJiP2v  -  PlvP2^i)ii  +  {P\vP2p  -  P\pP2v)^l 

+ {pipP2,.  -  pi^P2pyst]  ID ,  (S.4) 

where  pi  =  p-\-k  and  p2  =  p  —  k.  To  solve  the  model  equation  under  these  approxi- 
mations one  must  substitute  equations  {B.l)  to  {BA)  into  equations  (2.66)  through 
(2.69)  and  integrate.  This  can  be  accomplished  using  the  following  functional  rela- 
tion. Let  Xq  =  [x  ^  qf  —  {x'^  -\- 2x  ■  q -\-  q^).  Then  it  follows  that 

U^[x-'F{x^)\  =AF"  =G{xq)  iB.5) 

and  therefore, 

FM  =  I  ["  dx'^  £  rf<G«)  •  (5.6) 

Using  this  relation,  one  can  easily  show  that,  if 

nH(xJm^)  =  -^  ,  (5.7) 


Xg  —  m^ 
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then 
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$(x,/m2)  =  ^(l-^)ln(^-l).  (B.S) 

'1 


4  Xn  'TV? 


Also,  if  one  defines 

X^^{p  +  ikf,     ^^  =  (m'-P)  +  P^2     ^^^     Z^=X^/fi^  (5.9) 

then 


1  1  r       dc 


-ii: 


With  this  definition  it  follows  directly  that  if, 

1 


{pI  -m^)(pl  -m2) 


s 


It  is  easily  shown  that  one  can  write 


so  that,  if 


In 


-2  2 


P1-P2 

then 

1    /•'    ,..,„.       1    /•' 


then 


(5.10) 


then 


y'^^ln(Z^-l).  (S.12) 


(5.13) 


°'^6  -      .2    T2  '    '  (^-14) 


•^^  =  ^  /  ^  dmZ^)  =  IJ  ^  dCil  -  ^)  ln(Z^  -  1)  .  (5.15) 

Similarly  if 

□'*c  =  ^    ^im'-i,^)ln{Z^-l),  (5.16) 


<^c  =  ^  /',  ^^  {^^(^^)  "  (^  "  ^^^^(^^  "  ^)}  (^•^'^) 
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where  L2  =  -  J^  ln{l  -  z')^  is  the  dilogarithm.    The  dilogarithm  satisfies  the 
integral  relation 

/    L2{z')dz'  =  zL2{z)  +  (^  -  1)  ln(l  -z)-z  .  (5.18) 

Jo 

With  these  relations,  it  is  possible  to  solve  equations  (2.66)  to  (2.69)  as  follows: 
Equation  (2.66)  with  (B.l)  on  the  right-hand  side  becomes 

(p|  -  m^)(p^  -  m2)  dpx     pI-pI  ^         ' 

Integrating  Eq.(5.19)  using  Eq.(B.15)  gives 

Next,  substituting  Eq.(5.2)  into  Eq.(2.67)  yields 

^''''-'n'  (p1-m2)(pl-^2) -^1  {B.21) 

where  E^  satisfies  the  differential  equation 


<9p/.   [  (Pi  -  m2)(p2  _  ^2) 

One  can  easily  verify  that  the  following  relation  holds  true: 

(m^-pi-P2)^^  +  (pi^+P2^)p^  __\_d_  fpi^+p2M  j^  /^Pi  -  m2 


(S.22) 


(p?  -  m2)(p2  _  ^2)  2  apA  I  {p\  -  Pi)       \pI  -  m2 

Therefore,  C'^  can  be  written  in  the  form 
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where  C  and  C^,  satisfy  the  differential  equations 

□^C.  =  24i?lf^l"CI^)+?i}.  (B.25) 

\^{pI-P2)  \P2-^J         ^Pi^\ 


U-C  ^2e\,  S"^'  'J:-  ^'K.+  ^^  In  ^^^^^' 


(pI  -  m2)(p2  _  rn2)       pj-pl       \pI  -  m2 
1  1  4P 

=  eS  -735 ^  -  735 ^  + 


(pf-m2)       (p2_r7i2)       (pj  -  m'^){pl  -  m^) 
+  3^  in  (4^)},  (B.26) 

Pi-P5     Vp2-'^/J 


and 


(pf  -  m2)(p2  _  m2)       pj  -pl       \pl-m'^J' 
Using  Eqs.  (B.25)  to  (5.27)  it  is  found  that 

°'{'^-^)=''  '^•"'^ 

which  implies  that  C  =  -q-^.  This  also  imply  that 

d^C^,  =  0  .  (S.29) 

This  relation  is  valid  in  general  for  the  model  equation  as  can  be  seen  from  a  direct 
application  of  d'^  to  Eq.(2.67).  From  equations  {B.8)  and  (B.12),  one  obtains 


or 


il-X'-k-Sky^'^-'''^ 


c  =  -J  '''  ^^      ^       '^^ 


-(l-f)ln(§-l)-(l-f)..(§-lU.      (-30) 
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and  from  equation  (5.16)  with  Ct"  =  C  -  k^'Ccki'/P,  one  finds 


Pt^        In  f  ^1-"^    ]  +  ^ 


=  2e 


9p 


(5.31) 


This  follows  from  equations  (B. 12), (5. 15)  and  (5.22)  from  which 

I    [^   dE  f           m?         \ 
El  =  2m^$a  -  2$6  =  —  /     —     Z^  + 1     ln(Z^  -  1)  +  constant  .  (B.32) 

4  j-i  ^5  V         ^i       ) 

If  we  define  the  function  e  such  that  £?^  =  ep'^,  then  it  follows  that  e  —IdExjdX^^ 
and  thus 

/•I  r  1  /'™,2         ,,    \    r  1  \„(  r7  ^\\  ^ 

.(5.33) 
Using  Eq.(j5.16)  one  can  calculate 

4  =  -\^j^,  j\  ddm'  -  h){^2{Zs)  +  (1  -  ^)ln(Z(  -  1)} 

It  also  follows  from  Eqs.  (B.16)  and  (B.18)  that 


e  = 

"2/_/4m2£- 

(m^  -//^) 
-1)              /^i 

1 

In(Z^-l)' 
^1         J 

^2        _2 


D^PC,  =  2e  <!  -  In  (  ?^-^  ]^k-  VJSi 


4       \P2  ~  ^ 


i-^  —  m2 


(5.35) 


and 


□^(^%-ep-fcSi)  =  eln 


p\  —  m^ 


{pj  -m2)(p2  -^2) 


Therefore 


1 

4' 


^5.36) 


(5.37) 
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Using  Eq.(5.34)  with  (S.37)  gives 


'' = r^^  L  jky  -''A'-{'-^^  '^^''  -  '^] 


e 


1     f  p  •  k 


'^L 


'  di 


•(m2-;,^)ln(Z^-l) 


(5.38) 


4'v  ^'  y^^y-iA^f^e' 

Using  Eqs.  (5.24)  and  (5.38)  and  the  fact  that  the  function   v  multiphes  8^  in  the 
definition  of  C^,  we  finally  find  that 

It  also  follows  from  Eqs.(5.38)  and  (2.77)  that  if 


CkP^  =  k^'Ci  =  k  ■  (U2PI  +  UiP2)p^  , 


(5.40) 


then 


Cfc  =  -ep-  k 


■1  f^i^^ 


p^(Z^  -  1)       ^^Z^ 


HZ^  -  1)     .        (5.41) 


The  corresponding  equation  for  Cp  is 


Cpp^  =  p'^C^  ^[v+p-  {U2PI  +  UiP2)]  p^ 


(5.42) 


and  thus 


C    =-et      ^^     [(m^-/.g)(m^-P) 


1  .   ln(^^  -  1) 


-  fc^  ln(^5  -  1) 


(5.43) 
Using  Eqs.  (5.41)  and  (5.43)  it  is  possible  to  solve  for  the  functions  Ui  and  U2 
simultaneously  leading  to  the  expressions 

{Cp  -  v)pi  ■  k  -  Ckpi  ■  p 


Ui  = 


and 


U2 


2[p2p  _(p.fc)2] 
CkP2  -p-jCp-  v)p2  ■  k 

2[pH^  -ip-ky] 


(BAA) 


(5.45) 
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1  p-k 


or,  in  terms  of  the  (^-integral  transform, 

2 


ln(Z^  -  1) 


(m^  -  ;U^)-^  m 7T  +  2 


(5.46) 


and 


"^-rL]^y-^4'^[^r'h''-'^] 


1   p-k 
~  4^T2~ 


f      di         ,  f(Z^-2)      ,ln(Zg-i: 


^l^z.^ 


{z,  - 1] 


{BAl) 


Let  us  now  solve  for  the  functions  vi,V2  and  U3.    To  do  this,  we  substitute 
Eq.(B.3)  into  Eq.(2.68)  which  gives 


(pI  -rv?\pl  -m?) 


where  we  have  used  the  relation, 


□^C,\  =  elTDl  =  2meCSik  ■  □  -  M^);^;  _  „.)(g.  _  ^.)  ■  i^A,) 


From  Eqs.(5.12)  and  (5.49)  it  follows  that 


d'^C^^  =  2meC6ik  ■  □  -  ^,5^)  /    -^  ln(Z^  -  1)  (S.50) 

J-i  M«^J^ 


and  hence,  we  can  define  the  vector  functions 


1   -2.A  r  d( 


F^  =  Jf^a^C^^  =       ^mek'p 


liu 


1  /^«^«  L(^«  - 1)    ^^ 


'rn-i'"(^^-^= 


and 


f;=p''9'C  =  -i-*^pVl^ 


Ii^z-yA 


us^-i'"'"-^' 


(5.51) 


(5.52) 
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From  Eqs.(5.48)  and  (5.50),  we  find  that 
1 


D^C^^  =  -me 


{l^k  ■  U  -  K~d^)d,  -  {8lk  ■  U  -  k,~^^)^^   f  -\  \n{Z^  -  1) 


-  2me(kJ^  -  k,6^)— , 

and  thus,  using  Eq.(S.17),  one  obtains 


(5.53) 


C 


ixv 


+  jme{S^k,-S^k^) 


'     di 


H^i 


ln(Z^  -  1)  . 


(S.54) 


Notice  also  that 


{Slk-U-  k^d^)  j    di{L2{Z^)  +  (1  -  Z^)\n{Z^  -  1)}  = 


-/ 


1  H^i 


K^  ■  iP  +  m-P^k,\  <!  1  +  (^  -  l)ln(Z^  -  1) 


(5.55) 


Therefore,  the  function  u^  is  given  by 
Us  =  -me  /     — —r  ^ 


:i--)ln(Z^-l)-2 


(5.56) 


This  expression,  together  with  the  scalar  functions  fk  and  fp  allows  the  defini- 
tions, F^  -  fkp^  and  F^  =  fpp^,  in  Eqs.(5.51)  and  (5.52)  and  are  sufficient  to 
determined  the  functions  vi  and  V2.  The  calculation  is  straightforward  although 
somewhat  tedious.  The  final  form  of  the  functions  Vi  and  V2  expressed  in  terms 
of  the  (^-integral  transform  can  be  written  as  follows:  let  us  define  three  subsidiary 
functions  of  Z^  (  Ji,  I2,  /s) 


h 


l-(l--)ln(Z^-l) 


^2  =  -^  +  (-^--3)ln(Z^-l), 


(5.57) 

(5.58) 


and 


h  =  -Hz^-i). 

With  these  definitions,  the  functions  Vi  and  V2  are  given  by 
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(S.59) 


v^=-  -me  /     -^  [2/i  +  h] 
1-1  l^i 


8 


and 


V2 


+  -mel  4[p-^-(i-eK2-Pe(i-eK2]  (s.eo) 

=    \me  f  ^  [2/i  +  /3] 
o        7-1  ^^ 


+ 


-me         4[p.fc(l  +  e)J2  +  P^(l  +  eK2] 
^       J-l  /^^ 


(5.61) 


Finally,  the  equation  for  C^^^  is 


PltiP2u  -P\uP2^ 


(5.62) 


(pj  —  m?){p\  —  m?) 

{P\uP2p-pipP2u)Sl  +  {p\pP2^  -P\nP2p)Sl 

{p\  —  m'^)(p2  —  m^) 

Since  only  one  function  needs  to  be  evaluated,  one  can  try  to  solve  the  the  following 
4-vector  equation  without  any  loss  in  generality,  namely. 


^\p'iP2d''Ci,p)  =  -^^^d^ip'iP2Kp) 


[B.63) 


obtained  from  Eq.(5.62)  by  contracting  with  the  4- vectors  p^  and  P2-    Equation 
(5.63)  can  now  be  written  as 

(pIpI  -(Pi  ■P2f) 


np'iP2d''Ci.p)  =  -2e|a' 

+  8^ 


-2el  d 


{p1  —  m?)[p\  —  m'^) 

A  (Pi  •  P2)p2p  -  Plpip  +  (pi  •  P2)p\p  -  pIp2p 
(pl  -m'^Xpl  -m2) 

4A;V  1       ..  r  4PpV 


P 


(pj  -m2)(p2  _m2) 
16ePp^  -„~> 


d" 


{pi  —  m^)(pl  —  m?) 


/: 


(pj  —  m2)(p2  —  m^) 


A.^(Z^  -  1) 


(5.64) 
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Using  Eq.  (5.12)  one  has  that 

p^2^''c'^\p  =  -l^pp'  f  -^(1  -  hnz^  - 1)  •       (5.65) 

If  one  compares  this  equation  with  Eq.(5.20),  it  follows  that 

'^  '^  6  m 

The  scalar  function  ^3  is  given  by 

^'{vivici.,) = 5M^3  \m  -  (Pi  •  V2?\ } 

+  d"  {U3   [(Pl  ■  P2)(Plp  +  P2p)  -  (p?P2p  +  V\p\p)\  } 

=  -8Pv3p^  (5.67) 

as  terms  involving  derivatives  of  U3  disappear.    Thus  the  function  v^  is  related  to 

the  u-function  by  the  equation, 

1  u 
^3=0-'  (5.68) 

or  explicitly 


^3  =   i 


^^^/:ji{-^^}- 


For  reference  purposes,  all  functions  are  presented  in  Table  B-1.  Numerical  values 
for  the  functions  u,  ui,  U2,  vs  and  e  are  presented  in  Tables  B-2  to  B-6.  The  values 
of  the  variable  x  cover  the  perturbation  and  inner  asymptotic  regions. 
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Table  B-1    Perturbation  Solutions  to  the  Model  Equations 


Notation: 


X^  =  (p  +  fky,     fi^  =  (m^  -  P)  +  Fe  ,      Z^  =  X^/n^  , 


h  =  ^ 


1 


l-(l--)ln(Z^-i: 


,      /2  =  -^  +  (^-^)ln(Z^-l), 


'i 


zi     zi 


h  =  —\n{Z^-l) 


1         2)1- 


\  ln(Z^  -  \)di 


4^1  /_A^       Zi,       H^i) 


V  =  C 


e  = 


2  7-1      l/^H^?-i)  Mf 


1 ln(Zg  -  1) 


Z^(Ze  -  1) 


^1 


(5.70) 

(5.71) 
(5.72) 


Ui  =    -e 


1        /•'        <<«     ,      2 


1    p  •  /;; 


(m^-A««)|2  + 


L^5 


4^p  y-iA.pl ^"  "^^H(^r^    "" 


HZi  -  I] 


'^ 


U2 


-ri: 


im''-ii^)\2  + 


HZi  - 1)} 


1   p-k 


4^ 


[m 


-  ^h) 


{Zj  -  2)    ,  ^ln(Z^  -  1) 


+  2 


V'x   —  - 

4 


u  —  —em 
4 


4  "^   7-1  /^^^^  \^        z^        \ 

y_i  //^z^  1  z^      \  ■ 


(5.73) 


(5.74) 

(5.75) 
(5.76) 


Table  B-1  (Continued) 


76 


vi=--me         -^  [21,  +  h] 

+  1^'  f    §[P-  K^  -  Oh  -  Pe(l  -  0^2]   ,  (5.77) 

V2  =    \me  I    —  [2/1  +  73] 


+  7"^^/    4  [^•^'(1  +  0^2  +  ^^(1  +  0^2]   .  (5.78) 


=  i"^/.  ^  ['  -  i""'^^  -''-'}■  (^-^^ 
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Table  B-2    The  Perturbation  Functions  u,  ui,  i<2,  Vs  and  e  as  Functions 

2 

of  the  Variable  y  for  x  —  1  and  A;  =0.1 


_y_ 


tf 


(xlO^ 


ui(xlO^ 


U2(xl0^  )  v^jxlO^      )       e   (xlO°  ) 


1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

10.0 

20.0 

30.0 

40.0 

50.0 

60.0 

70.0 

80.0 

90.0 

100.0 

200.0 

300.0 

400.0 

500.0 

600.0 

700.0 

800.0 

900.0 

1000.0 


-3.0219 

-1.9098 

-1.3252 

-0.9509 

-0.6881 

-0.4933 

-0.3436 

-0.2254 

-0.1304 

-0.0527 

0.2891 

0.3686 

0.3847 

0.3802 

0.3684 

0.3543 

0.3398 

0.3256 

0.2270 

0.1755 

0.1430 

0.1200 

0.1027 

0.0891 

0.0782 

0.0693 

0.0620 

0.0276 

0.0165 

0.0112 

0.0082 

0.0064 

0.0051 

0.0042 

0.0035 

0.0030 


-5.4450 

-5.4118 

-5.1836 

-4.9085 

-4.6286 

-4.3587 

-4.1044 

-3.8671 

-3.6469 

-3.4428 

-2.0523 

-1.3236 

-0.8929 

-0.6161 

-0.4275 

-0.2934 

-0.1948 

-0.1205 

0.1363 

0.1795 

0.1907 

0.1888 

0.1789 

0.1648 

0.1489 

0.1331 

0.1182 

0.0313 

0.0024 

-0.0093 

-0.0145 

-0.0171 

-0.0182 

-0.0187 

-0.0188 

-0.0186 


5.1248 

5.5385 

5.5005 

5.3192 

5.0860 

4.8372 

4.5889 

4.3486 

4.1197 

3.9037 

2.3685 

1.5404 

1.0481 

0.7318 

0.5166 

0.3637 

0.2515 

0.1670 

-0.1252 

-0.1752 

-0.1885 

-0.1877 

-0.1783 

-0.1644 

-0.1488 

-0.1331 

-0.1183 

-0.0315 

-0.0025 

0.0091 

0.0144 

0.0170 

0.0182 

0.0186 

0.0187 

0.0186 


-1.0073 

-0.6366 

-0.4417 

-0.3170 

-0.2294 

-0.1644 

-0.1145 

-0.0751 

-0.0435 

-0.0176 

0.0964 

0.1229 

0.1282 

0.1267 

0.1228 

0.1181 

0.1133 

0.1085 

0.0757 

0.0585 

0.0477 

0.0400 

0.0342 

0.0297 

0.0261 

0.0231 

0.0207 

0.0092 

0.0055 

0.0037 

0.0027 

0.0021 

0.0017 

0.0014 

0.0012 

0.0010 


-0.6688 
-0.6167 
-0.5797 
-0.5500 
-0.5249 
-0.5031 
-0.4838 
-0.4665 
-0.4508 
-0.4365 
-0.3387 
-0.2821 
-0.2441 
-0.2164 
-0.1951 
-0.1781 
-0.1643 
-0.1527 
-0.0927 
-0.0684 
-0.0548 
-0.0460 
-0.0399 
-0.0353 
-0.0317 
-0.0288 
-0.0265 
-0.0150 
-0.0107 
-0.0084 
-0.0069 
-0.0059 
-0.0052 
-0.0046 
-0.0042 
-0.0038 
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Table  B-3   The  Perturbation  Functions 


U,    Ui,    U2, 
2 

of  the  Variable  y  for  a;  =  10  and  k    =0.1 


U3  and  e  as  Functions 


X 


(xlO^ 


uiJxlO^ 


U2(xl0^  )   t;3(xlO^  )   e  (xlO"  ) 


1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

10.0 

20.0 

30.0 

40.0 

50.0 

60.0 

70.0 

80.0 

90.0 

100.0 

200.0 

300.0 

400.0 

500.0 

600.0 

700.0 

800.0 

900.0 

1000.0 


4.1158 
3.9866 
3.8658 
3.7525 
3.6458 
3.5452 
3.4500 
3.3597 
3.2739 
3.1922 
2.5381 
2.0595 
1.6313 
1.2907 
1.1057 
0.9760 
0.8780 
0.8004 
0.4456 
0.3175 
0.2482 
0.2036 
0.1721 
0.1484 
0.1301 
0.1155 
0.1036 
0.0489 
0.0307 
0.0218 
0.0165 
0.0130 
0.0105 
0.0086 
0.0072 
0.0061 


2.2057 

-2.1843 

2.9492 

-2.9287 

3.4246 

-3.4049 

3.7571 

-3.7383 

3.9982 

-3.9802 

4.1754 

-4.1582 

4.3058 

-4.2893 

4.4006 

-4.3848 

4.4675 

-4.4524 

4.5123 

-4.4979 

4.3514 

-4.3420 

3.7814 

-3.7753 

2.9726 

-2.9689 

2.2549 

-2.2533 

1.9353 

-1.9349 

1.7255 

-1.7259 

1.5713 

-1.5722 

1.4506 

-1.4518 

0.8852 

-0.8864 

0.6699 

-0.6708 

0.5547 

-0.5553 

0.4782 

-0.4786 

0.4193 

-0.4196 

0.3705 

-0.3708 

0.3290 

-0.3293 

0.2935 

-0.2937 

0.2629 

-0.2631 

0.1068 

-0.1068 

0.0535 

-0.0536 

0.0284 

-0.0284 

0.0141 

-0.0141 

0.0051 

-0.0051 

0.0011 

0.0011 

0.0055 

0.0055 

0.0087 

0.0087 

0.0112 

0.0112 

1.3719 

-0.1073 

1.3289 

-0.1045 

1.2886 

-0.1020 

1.2508 

-0.0996 

1.2153 

-0.0974 

1.1817 

-0.0953 

1.1500 

-0.0934 

1.1199 

-0.0916 

1.0913 

-0.0899 

1.0641 

-0.0883 

0.8460 

-0.0772 

0.6865 

-0.0730 

0.5438 

-0.0924 

0.4302 

-0.0557 

0.3686 

-0.0447 

0.3253 

-0.0381 

0.2927 

-0.0335 

0.2668 

-0.0301 

0.1485 

-0.0157 

0.1058 

-0.0110 

0.0827 

-0.0085 

0.0679 

-0.0070 

0.0574 

-0.0060 

0.0495 

-0.0052 

0.0434 

-0.0046 

0.0385 

-0.0042 

0.0345 

-0.0038 

0.0163 

-0.0021 

0.0102 

-0.0015 

0.0073 

-0.0011 

0.0055 

-0.0009 

0.0043 

-0.0008 

0.0035 

-0.0007 

0.0029 

-0.0006 

0.0024 

-0.0005 

0.0020 

-0.0005 
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Table  B-4   The  Perturbation  Functions  u,  ui,  1/2,  v^  and  e  as  Functions 


2 

of  the  Variable  y  for  x  —  10^  and  k    —  0.1 


u 


(xlO*^     )      ui(xlO^     )      U2(xl0^     )      vzJxlO^     )      e  (xlO^     ) 


1.1 

3.3953 

0.1657 

-0.1665 

1.1318 

1.2 

3.2917 

0.2229 

-0.2235 

1.0972 

1.3 

3.1954 

0.2600 

-0.2606 

1.0651 

1.4 

3.1055 

0.2865 

-0.2871 

1.0352 

1.5 

3.0215 

0.3062 

-0.3068 

1.0072 

1.6 

2.9428 

0.3212 

-0.3217 

0.9809 

1.7 

2.8687 

0.3328 

-0.3332 

0.9562 

1.8 

2.7990 

0.3417 

-0.3421 

0.9330 

1.9 

2.7331 

0.3486 

-0.3490 

0.9110 

2.0 

2.6709 

0.35.38 

-0.3542 

0.8903 

3.0 

2.1921 

0.3625 

-0.3628 

0.7307 

4.0 

1.8762 

0.3448 

-0.3450 

0.6254 

5.0 

1.6496 

0.3233 

-0.3234 

0.5499 

6.0 

1.4780 

0.3027 

-0.3028 

0.4927 

7.0 

1.3427 

0.2841 

-0.2842 

0.4476 

8.0 

1.2331 

0.2676 

-0.2676 

0.4110 

9.0 

1.1420 

0.2529 

-0.2529 

0.3807 

10.0 

1.0651 

0.2398 

-0.2399 

0.3550 

20.0 

0.6579 

0.1614 

-0.1614 

0.2193 

30.0 

0.4882 

0.1248 

-0.1248 

0.1627 

40.0 

0.3913 

0.1034 

-0.1034 

0.1304 

50.0 

0.3269 

0.0889 

-0.0890 

0.1090 

60.0 

0.2805 

0.0781 

-0.0782 

0.0935 

70.0 

0.2452 

0.0696 

-0.0696 

0.0817 

80.0 

0.2176 

0.0627 

-0.0627 

0.0725 

90.0 

0.1953 

0.0569 

-0.0569 

0.0651 

100.0 

0.1770 

0.0520 

-0.0520 

0.0590 

200.0 

0.0904 

0.0273 

-0.0273 

0.0301 

300.0 

0.0603 

0.0183 

-0.0183 

0.0201 

400.0 

0.0452 

0.0137 

-0.0137 

0.0151 

500.0 

0.0361 

0.0110 

-0.0110 

0.0120 

600.0 

-0.1990 

-0.0179 

0.0179 

-0.0663 

700.0 

-0.2031 

-0.0381 

0.0381 

-0.0677 

800.0 

-0.2057 

-0.0536 

0.0536 

-0.0686 

900.0 

-0.2071 

-0.0672 

0.0672 

-0.0690 

1000.0 

-0.1942 

-21.3921 

21.3921 

-0.0647 

-9.6858 
-9.3891 
-9.1132 
-8.8560 
-8.6154 
-8.3899 
-8.1779 
-7.9782 
-7.7897 
-7.6114 
-6.2410 
-5.3369 
-4.6888 
-4.1979 
-3.8114 
-3.4979 
-3.2379 
-3.0182 
-1.8564 
-1.3738 
-1.0994 
-0.9179 
-0.7874 
-0.6886 
-0.6112 
-0.5490 
-0.4980 
-0.2556 
-0.1712 
-0.1286 
-0.1029 
-0.1316 
-0.1100 
-0.0954 
-0.0845 
-0.0760 
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Table  B-5    The  Perturbation  Functions  u,  ui,  1x2,  ^3  and  e  as  Functions 


—2 


of  the  Variable  y  for  x  =  10      and  k    =  0.1 


u   (xlO^'^  )   ui(xlO^^  )   U2(xl0^^  )   t;3(xl0^^  )  e   (xlO^i 

-9.6795 
-9.3825 
-9.1064 
-8.8489 
I  1.5      2.9993      3.0253      -3.0253      0.9998      -8.6082 

I  1.6      2.9206      3.1727      -3.1727      0.9735      -8.3825 

-8.1703 
-7.9705 
-7.7818 
-7.6035 
-6.2323 
-5.3278 
-4.6794 
-4.1884 
-3.8017 
-3.4882 
-3.2281 
-3.0083 
-1.8463 
-1.3639 
-1.0902 
-0.9098 
-0.7804 
-0.6827 
-0.6062 
-0.5448 
-0.4943 
-0.2545 
-0.1707 
-0.1283 
-0.1028 
-0.0857 
-0.0735 
-0.0643 
-0.0572 
-0.0514 


1.1 

3.3725 

1.6408 

-1.6408 

1.1242 

1.2 

3.2691 

2.2043 

-2.2043 

1.0897 

1.3 

3.1729 

2.5703 

-2.5703 

1.0576 

1.4 

3.0832 

2.8312 

-2.8312 

1.0277 

1.5 

2.9993 

3.0253 

-3.0253 

0.9998 

1.6 

2.9206 

3.1727 

-3.1727 

0.9735 

1.7 

2.8467 

3.2859 

-3.2859 

0.9489 

1.8 

2.7771 

3.3731 

-3.3731 

0.9257 

1.9 

2.7114 

3.4403 

-3.4403 

0.9038 

2.0 

2.6492 

3.4914 

-3.4914 

0.8831 

3.0 

2.1715 

3.5711 

-3.5711 

0.7238 

4.0 

1.8563 

3.3919 

-3.3919 

0.6188 

5.0 

1.6304 

3.1761 

-3.1761 

0.5435 

6.0 

1.4593 

2.9703 

-2.9703 

0.4864 

7.0 

1.3246 

2.7848 

-2.7848 

0.4415 

8.0 

1.2154 

2.6200 

-2.6200 

0.4051 

9.0 

1.1247 

2.4740 

-2.4740 

0.3749 

10.0 

1.0482 

2.3442 

-2.3442 

0.3494 

20.0 

0.6433 

1.5661 

-1.5661 

0.2144 

30.0 

0.4752 

1.2033 

-1.2033 

0.1584 

40.0 

0.3798 

0.9914 

-0.9914 

0.1266 

50.0 

0.3170 

0.8491 

-0.8491 

0.1057 

60.0 

0.2719 

0.7442 

-0.7442 

0.0906 

70.0 

0.2379 

0.6622 

-0.6622 

0.0793 

80.0 

0.2112 

0.5961 

-0.5961 

0.0704 

90.0 

0.1898 

0.5413 

-0.5413 

0.0633 

100.0 

0.1722 

0.4954 

-0.4954 

0.0574 

200.0 

0.0887 

0.2634 

-0.2634 

0.0296 

300.0 

0.0595 

0.1779 

-0.1779 

0.0198 

400.0 

0.0447 

0.1341 

-0.1341 

0.0149 

500.0 

0.0358 

0.1075 

-0.1075 

0.0119 

600.0 

0.0299 

0.0897 

-0.0897 

0.0100 

700.0 

0.0256 

0.0769 

-0.0769 

0.0085 

800.0 

0.0224 

0.0674 

-0.0674 

0.0075 

900.0 

0.0199 

0.0599 

-0.0599 

0.0066 

1000.0 

0.0179 

0.0539 

-0.0539 

0.0060 

81 


Table  B-6    The  Perturbation  Functions  u,  Ui,  ^2,  ^3  and  e  as  Functions 
of  the  Variable  y  for  x  =  10'^°  and  k    =  0.1 


y  u   (xlO^^  )   ui(xl0^4  )   tz2(xlO^^  )   ^3(xl0^^  )   e  (xlO^^  ) 

1.1  3.3725  1.6408  -1.6408  1.1242  -9.6795 

1.2  3.2691  2.2043  -2.2043  1.0897  -9.3825 

1.3  3.1729  2.5703  -2.5703  1.0576  -9.1064 

1.4  3.0832  2.S312  -2.8312  1.0277  -8.8489 

1.5  2.9993  3.0253  -3.0253  0.9998  -8.6082 

1.6  2.9206  3.1727  -3.1727  0.9735  -8.3825 

1.7  2.8467  3.2859  -3.2859  0.9489  -8.1703 

1.8  2.7771  3.3731  -3.3731  0.9257  -7.9705 

1.9  2.7114  3.4403  -3.4403  0.9038  -7.7818 
2.0  2.6492  3.4914  -3.4914  0.8831  -7.6035 
3.0  2.1715  3.5711  -3.5711  0.7238  -6.2323 
4.0  1.8563  3.3919  -3.3919  0.6188  -5.3278 
5.0  1.6304  3.1761  -3.1761  0.5435  -4.6794 
6.0  1.4593  2.9703  -2.9703  0.4864  -4.1884 
7.0  1.3246  2.7848  -2.7848  0.4415  -3.8017 
8.0  1.2154  2.6200  -2.6200  0.4051  -3.4882 
9.0  1.1247  2.4740  -2.4740  0.3749  -3.2281 

10.0  1.0482  2.3442  -2.3442  0.3494  -3.0083 

20.0  0.6433  1.5657  -1.5657  0.2144  -1.8464 

30.0  0.4758  1.1994  -1.1994  0.1586  -1.3654 

40.0  0.3818  0.9826  -0.9826  0.1273  -1.0958 

50.0  0.3210  0.8378  -0.8378  0.1070  -0.9212 

60.0  0.2780  0.7334  -0.7334  0.0927  -0.7980 

70.0  0.2460  0.6541  -0.6541  0.0820  -0.7060 

80.0  0.2211  0.5917  -0.5917  0.0737  -0.6344 

90.0  0.2010  0.5412  -0.5412  0.0670  -0.5770 

100.0  0.1846  0.4993  -0.4993  0.0615  -0.5299 

200.0  0.1044  0.2899  -0.2899  0.0348  -0.2996 

300.0  0.0743  0.2089  -0.2089  0.0248  -0.2132 

400.0  0.0582  0.1651  -0.1651  0.0194  -0.1671 

500.0  0.0481  0.1372  -0.1372  0.0160  -0.1382 

600.0  0.0412  0.1179  -0.1179  0.0137  -0.1182 

700.0  0.0361  0.1036  -0.1036  0.0120  -0.1035 

800.0  0.0321  0.0926  -0.0926  0.0107  -0.0922 

900.0  0.0290  0.0838  -0.0838  0.0097  -0.0833 

1000.0  0.0265  0.0767  -0.0767  0.0088  -0.0760 


APPENDIX  C 
ASYMPTOTIC  FORMS  OF  GREEN'S  PERTURBATION  SOLUTIONS 

In  order  to  obtain  the  asymptotic  forms  of  the  perturbation  solutions  derived 
in  Appendix  B,  we  must  observe  that  the  x-dependency  of  the  i^-integral  transform 
occurs  in  the  function  X^  —  p^  +  2p  ■  k  +  k'^ .  In  the  asymptotic  region, 


p^  ^  xy;  p-  k  —  \x\/ip'  -  1   .  (C.l) 

Consequently,  X^  can  be  approximated  by  the  expression 

where 

©(0  =  y  +  e\/y'-i  •  (C.3) 

Notice  also  that  Z^  —  1  can  be  written  in  the  form 

^e  - 1  =  (n)  #^  .  (C.4) 

where 


k^j  e+x' 


-  1  ,  (C.5) 


and  thus 


ln(Z^  _  1)  ^  In  ^  +  In  0(0  -  He  +  x')  •  (C.6) 

To  obtain  a  good  approximation  of  the  i^-integrals,  in  the  asymptotic  re- 
gion, it  is  necessary  to  keep  only  those  terms  in  the  integrand  that  have  the 
smallest  power  in  x~^ .    We  must  be  careful,  not  to  neglect  terms  such  as  Z7^ 
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against  Z^^lnZ^  due  to  the  slow  variation  of  the  logarithm.  To  avoid  confu- 
sion between  the  perturbation  solutions  in  the  asymptotic  region  and  the  actual 
asymptotic  solutions,  we  will  use  the  superscript  "pert"  when  we  refer  to  the  former 
solutions. 

Using  Eq.(S.20),  (See  also  Integral  Table  at  the  end  of  the  chapter)  it  follows 
that 

This  expression  is  not  accurate  enough  in  the  asymptotic  region  due  to  the  fact 
that  in  obtaining  the  perturbation  solutions  we  set  Ai  ~  A2  ~  m.  However,  by 
Eq.(3.61),  Ai  ~  v42  ~  x^'^ .  It  has  been  shown^^'^^  that  an  extended  solution  can 
be  constructed  if  we  replace  the  mass  m  by  the  average  of  Ai  and  A2 .  Assuming 
this  to  be  valid,  then 

uiP-^)  =  liA,+A,)e—^.  iC.8) 

The  form  of  the  function  u^^^"^')  in  the  asymptotic  region  can  be  obtained  from 
Eqs.(S.30)  and  (5.39).  From  Eq.(5.30),  we  see  that  it  is  possible  to  approximate 
the  function  C  by  the  form 

1      /•\    .„       ..,.      1 


m4)+m4) 


((7.9) 


and  therefore, 

v(P^^'^r.l+C^l  +  le\l-2xt^n-\l/x)  +  -jt^]   .  (C.IO) 

-    [  vy^  -  1  J 

Notice  that  the  largest  term  in  u(p<='"*)  is  x-independent;  all  other  terms  being  of 
order  0(l/x).  From  Eq.(B.33)  it  also  follows  that  the  function  e^^'^^')  is  given  by 


e 


1  /'  ^ 

2  7-1  X. 


(pert)  ^      1^     dC  InZ 


(C.ll) 
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The  form  of  the  functions  u)^^^  and  ^3  '^'"  can  be  extracted  from  Eqs.(5.46) 
and  (jB.47).  It  is  clear  that  the  largest  contribution  to  these  functions  comes  from 
the  integrals  being  multiplied  by  p  ■  ^.  From  these  equations,  it  follows  that  there 
is  a  simple  relation  between  the  two  functions,  i.e. 


{pert)  (pert)  J-        /  "C      /      2  \ 


=     —Ur, 


This  integral  is  elementary  (see  Table  C-1),  and  hence 


iC.12) 


u 


{pert) 


=  —u 


(pert) 


y  In  Z  1 

iy'  -  1)  "  7^^^ 


(CIS) 


We  now  turn  to  the  problem  of  finding  the  form  of  the  functions,  u^^*'"*^  ^  v^'^^' 
and  u/^^  ,  in  the  asymptotic  region.  These  functions  are  particularly  difficult  to 
calculate  since  we  must  split  the  terms  containing  the  function  ln(Z^  —  1)  as  pre- 
scribed in  Eq.(C.6).  Also  a  test  of  these  functions  by  substitution  in  their  cor- 
responding differential  equations,  show  that  there  is  not  a  simple  assumption  we 
can  make  regarding  the  transition  of  the  mass  m  from  its  perturbation  value  to  its 
asymptotic  form. 

The  function  u/^^  given  in  Eq.(.B.56)  has  the  following  form  in  the  asymp- 
totic region. 


u 


(pert) 


^ e 


'      dC     u 


^    -l-if^lZl 


2  +  ln(Z^-l)]   . 


(C.14) 


Using  Eq.(C.6)  and  the  integrals  found  in  Table  C-1  it  is  possible  to  write  U3  ''*^ 
in  a  somewhat  complex  form.  To  this  end,  let  us  define  the  functions  tj  and  T2  by 
the  equations 


ylnz 


and     r2  = 


P 


2(y2_l)  +  ^^2 


m 


Using  these  definitions,  Eq.(Cl4)  integrates  into  the  form, 


(C.15) 


(pert) 


-t  [1  -  \n{xlm')  +  n  +  2r2[x'tan-i(l/x)  -  n]]  jx' 


(ai6) 
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where  x  is  defined  in  Eq.(C.5). 

Tlie  form  of  vY^^     in  the  asymptotic  region  can  be  found  in  an  analogous 
manner  from  Eq.(B.60).  In  this  region  we  have 


(pert) 


r"/,{-;4"-'"'^^"'" 


2p-k 

,,2  72 


(l-0[2-ln(Z^-l)] 


Using  Eq.(C6)  and  defining  the  integrals, 


h  =  !]''({' Y  +  T  N-^/t')  +  ine{o]  - 


"1^(1-0 


2n  ■  k 

+  4^(i-0[inWP)  +  ine(e)] 


and 


(C.17) 


(CIS) 


it  is  possible  to  write  Eq.(Cl7)  in  the  form, 

vi'^"-'^  -  Imeih  -  h)  .  {C.20) 

Using  Table  C-1,  it  is  found  that  the  integral  Ii  can  be  written  in  the  form 


/i  =  [In  -  2p  ■  fc(Ii3  -  /14)]  [ln(a;/m2)  -2]+2p-  kl^s  ,  {C.21) 


where 


L 


11 


hi 


'  Cd^        2 


7-1  ^1 


Inz 


(y^  - 1)    V^f^ 


(C.22) 


Je         2 


,      /i4  =    /     1;^  =  7^  ,      (C.23) 


X|       x^ 


and 


tl5 


/: 


^   (l-e)ln0(eK 


^1 


'Vy^^ 


^-(-^  +  \/y^-l)lnJ    .        (C.24) 
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Substituting  Eqs.(C.22)  to  (C.24)  into  Eq.(a21)  gives 


2z  Hx/¥) 


(C.25) 


The  integral  I2  given  in  Eq.(Cl9)  can  be  split  into  three  integrals  of  the  form, 


1  r 


hi  +  I22  Vy^  - 1  ~  hsyy^-l 


{C.26) 


where 


I21 
-^22 


and 


^23  — 


J-i       0(0        ' 

^  f  He+x')dc 
'  cHe  +  x')d^ 


-1 


02(0 


(C.27) 
(C.28) 

(C.29) 


The  integral  in  Eq.(C.27)  cannot  be  expressed  in  terms  of  a  finite  number  of  ele- 
mentary functions.  Fortunately,  it  will  not  be  necessary  to  evaluate  this  integral  in 
order  to  define  v^'^'^  (or  ^2  ^"^  )  since  the  integral  I23  can  be  expressed  in  terms 
of  the  integral  I21.  To  show  this,  we  must  integrate  by  parts  Eq.(C29),  i.e. 


hs  =  (L. 


22 


e=-i 


-/ 


hidC  . 


(C.30) 


In  general,  we  have 


GW  =  / 


\n(t^  +  a^)dt 
(ct  +  df 


I 


tdt 


lln(^2  +  a2)       2 

c    {ct  +  df    ^cj   (t^  +  a^)(ct  +  d)  ' 


(C.31) 


anc 


/ 


tdt 


{f"  +  a^){ct  +  d)       a^c^  +  d^^2 


{-aln(t^  +a^)  -  dln{ct  +  d)  +  actan~^{t/a)}  . 


(C.32) 
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Setting  a  =  X,  c  =  ^xp'  -  1  and  <i  =  y  in  Eq.(C.31)  and  evaluating  the  integral 
gives 

/22  =  21n(l  +  x')  +  4r2[xtan-i(l/x)  -  t^\  (CSS) 

where  Tj  and  T2  are  defined  in  Eq.(C.14).  To  complete  the  evaluation  of  the  integral 
I22,  we  must  integrate  once  again  Eq.(C.31).  This  integration  is  elementary,  and 
yields  the  result, 

J-\       Jo 


'^      ln{t^  +  a'^)dt  _  /21 


(cf  +  dy         Vy^^ 


+     f^  0^^^  +  ^'')  +  2xtan-^(l/x)  -  2x2}    .        (CM) 

Vy  -1 


Multiplying  Eq.(C31)  by  t  and  evaluating  the  expression  at  t  =  ±1  gives 

2 
-2^r2yVy'-lln(mVP)  •  (^.35) 


tG(i) 


t=^+i  2 


t=-i  ''^ 


Using  these  results,  we  find  that  the  integral  I2  can  be  written  in  the  form, 

h  -  ^^'~  ^  [I22  -  r^{y)]  ,  (CM) 

where  the  function  T3  is  defined  by 

"■^(y)  =  -^^P=f  {ln(mVP)/r2  +2xtan-i(l/x)  -  2ti  }   .  {C.37) 

Notice  that  the  integral  /21  has  disappeared. 

Finally,  substituting  Eqs.(C.36)  and  (0.25)  into  Eq.(C.20)  yields  the  desired 
result  for  the  function  vj^^^  ' ,  namely 

where 

H{x,y)  =  l-\n{x/m^)  +  ^^L^  +  2T2[xtSin-\l/x)-r,]  .  (0.39) 

vy  -1 

A  similar  procedure  can  be  used  to  show  that  the  function  V2  is  given  by  the 

expression. 


V. 


(pert)  -I-     —  /  ■.  ^ 


2 


=  +-6i?(x,y) .  (C.40) 


4       '   ""'    X 


For  convenience,  all  asymptotic  forms  are  presented  in  Table  C-2. 


88 


Table  C-1    Table  of  Useful  Integrals. 


Notations: 


I 3;  u 


x^/y'^  -  1 


m'- 


S^ 


(xz  —  m^) 
(xz~^  —  m^ 


;      R  =  {xz  —  m?){xz    ^—m^)/m^;      x 


1; 


T2  = 


m2(y2_l)  +  p    ■ 


1. 


dC 


2lnz 


'■■I 


/- 


a  +  K)       xyjy'^  -  1 
di  2 


a +  60       a;2  ' 
idi  2 


a  +  ^0        xy?/2  _  1 


edi 


-2y 


a  +  hi)        x{t  -  1) 
de  2 


1- 


ylnz 
y  Inz 


a  +  6^2       ^2 

e^e       2 


edi 


\nz 


a  +  bi)       x2(j/2_i) 
di  2y 


(y^  - 1)     ^^^^ 


[l  +  y2-yln. 


a  +  6^)3       x^ 
di 


10.  /    ln(a  +  ^0^^  =  -2  +      ^j^^"^      +  2\n{x/m^) 
J-i  a^vy^  -  1 


a  +  bi)(a  +  I3i)        x^        m^x^y^  -  1 


[In5-21n2] 
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Table  C-1  (Continued) 

'  ln(a  +  /30dC        1  [ln(xz/m2)]2  -  [ln{xz-^ /m^)f 


10 


11 


12 


■/; 


(a  +  60 


x\^y'^  -  1 


r  ln(a 


+  Z'0 

^  ln(a  +  60(^C 


0  . 


J-i     {a 


'  (  Ha  +  bOdC 


1- 


ylnz 


13.  /'  iif 

J-1  (« 


2y 


+  60' 


+  21n. 


14.  /'  if 

J-1      [a 


+  bO' 


21n(l  +  x')  +  4r2 


1  + 


Xtan-I(l/X)- 


y  Inz 


V^ 


Table  C-2   Asymptotic  Forms  of  the  Perturbation  Solutions 
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Notation: 


y\nz 


P 


2yr2 


m?{y^  -  1)  +  P   ' 
^3^^)  = ^1=^  {ln(mVP)/r2  +  2xtan-^(l/x)  -  2ti}    ; 


H{x,  y)  =  1  -  ln(x/m^)  + 


y  ln5r 


2r2[ri-xtan-i(l/x)]  ;      x  = 


'm' 


1 
4' 


^iP-rt)  ^i+c^i^-;2-  4xtan-^(l/x)  + 


2ylnz 


,ipert) 


(pert)  ,    J- 


(pert)  J- 

^  2 


(peri) 


ln^ 


1  In^ 


^3  ~  -e 


(pert) 


2    x-v/y^  _  j^ 
3  ,  .  .    ,  In  z 


u^p--^>  r^  :-(A,  +  A2)e 


Vy'-i  ' 


.?-'U+iei7(x,y)i. 


P       "^  ■ 


(C.41) 

(C.42) 

(C.43) 

(a44) 

(C.45) 

(C.46) 

(C.47) 
(a48) 

(a49) 


APPENDIX  D 
ASYMPTOTIC  BOUNDARY  CONDITIONS 

To  obtain  the  boundary  conditions  to  be  used  with  the  differential  equa- 
tions derived  from  the  model  equation,  we  use  a  technique  developed  and  used 
by  A.  Broyles^°  in  connection  with  the  determination  of  the  electron  propagator 
functions.  This  technique  is  equivalent  to  the  use  of  Green's  theorem  by  which  it  is 
possible  to  convert  certain  volume  integrals  into  surface  integrals. 

In  order  to  find  the  boundary  conditions  for  the  vertex  function,  we  start  with 
the  definition  of  the  vertex  integral  equation  as  given  in  Eq.(2.34),  namely, 

L\puP2)^xt  +  ^^l^jD,Ap-q)rE.^%quq2)d'q  (D.l) 


where 


and 


D^„{k)  =  Zz  [-g^^  +  k^k^k-^]  k-^  .  {DA) 


We  have  also  defined  the  4- vectors  pi  =  p  +  k  ,  p2  =  p  —  k  which  represent  the 
outgoing  and  incoming  electron  momenta  respectively.   Substituting  Eq.(D.4)  into 

Eq.(D.l)  gives 

r  {PUP2)  -1+j^J  {-(-— )i  + ^3-)i }d  ^  ■      iD-5) 
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Using  the  identity, 
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P^Pi^P  "^  =  --d^d^lnp'^  +  -g^^p   2  ,  (£).6) 


the  last  term  in  Eq.(£).5)  can  be  written  as 


rt^r^^^^^-^im-^   <-) 


where 


^ 


Substituting  Eq.{D.7)  into  (D.5)  yields 

We  can  apply  )^  to  Eq.(D.9)  with  the  result, 

We  would  like  to  expresss  the  integral  Q^  so  that  the  p-dependency  appears  only  in 
the  form  (p  —  q)~'^ .  Notice  that  from  Eq.(£).8), 

^5=-|[Dja,^ln(p--g-)2]Z^^A 

+  J{niHp-qf]d,X''d'q,  (DAI) 

where  D^  stands  for  the  D'Alambertian  operator  in  the  momentum  q  (see  also 
Eq.(2.35))  and  dq'^  —  d/dq^.  In  general,  a  subscript  q  will  be  used  to  denote 
derivatives  with  respect  to  q.  The  first  integral  in  Eq.(D.ll)  can  easily  be  obtained 
and  is  given  by 

Jdq^[[U'^Hp-qf]F'-^]d\=<f{]Pllnip-qf]F'-^]da^  {D.12) 
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where  da^,  is  the  outward  drawn  normal  to  the  volume  d'^q.  The  identity, 

nHnp^  =  4p'\  (D.13) 

can  be  used  to  show  that 

Also,  from  the  definition  of  )^,  it  follows  that 

=  -2'"2./p,.(p-«n£"'A. 

It  is  possible  to  transform  Eq.(D.15)  using  the  identity, 

d,Jp  -  qy'F^"]  =  [d,,{P  -  q)-']F"^  +{p-  q)-'d,^F^^  .  (D.16) 

Using  Eq.(£).16)  in  Eq.(D.15)  yields  the  result 

Substituting  Eq.(D.17)  and  (D.U)  into  Eq.(D.lO)  we  obtain 

1     ie     (fd^-fF^"        rf-f  F^^d^q^ 


or,  grouping  all  integrals  over  the  volume  J^g, 


(D.15) 


(D.18) 


+  (27r)4v    (p-gT    ^J{p-q? 
We  can  reduce  Eq.(£).18)  to  a  differential  equation  by  using  the  identity 


i^     UI^^1..E-''       IF^'"^.'  ^^-''^ 
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□V'  =  i{2wY8\p)  .  {D.20) 

If  we  exclude  an  infinitesimal  volume  around  the  point  p  =  q,  there  will  be  no 
contributions  from  the  surface  integrals  in  Eq.(D.19).  Hence,  if  we  apply  the 
D'Alambertian  D^  to  Eq.(£>.19),  we  obtain 

□')^r\pi,p2)  =  Kd.F^"  +  Ifx^F^n  .  (D.21) 

Notice  that  this  is  the  same  equation  obtained  in  Chapter  II  by  an  independent 
method  (see  Eq.(2.50)). 

In  order  to  obtain  the  boundary  conditions  to  the  differential  equation  given 
in  Eq.(£).21),  we  integrate  Eq.(£).8)  by  parts  excluding  the  point  (p  =  q)  to  obtain 

=  -f  j>\n{p  -  qfjt'da,  +  IJ   f  \n{p  -  qfdg^F^'dS  .  {D.22) 

We  can  rearrange  Eq.(D.21)  to  obtain  the  form, 

ed.F^"  =  D^^f^r^  -  lefj^F^"  .  (D.23) 

Substituting  Eq.(Z).23)  into  ^(^.{D.22)  gives 

e^^-ef  j  \n{p  -  qfF^^'da,  +f  f  ln(p  -  qfUp^T^^d^ 

+  \elJ  fHp-Qf^g%E.^"d\  ■  (DM) 

Integrating  by  parts  the  last  two  terms  in  Eq.(I?.24)  gives 

eQ  =  -ef  I  ln(p  -  qfF^^da,  +  f<j>  ln(p  -  qfU]d^V^ 

-  iJ   UiJ,  ln(p  -  qf\U'^Y}d\  -  \elJ  /ln(p  -  qfd^j^^F^'' 

+  \e1/  j[f,  ln(p  -  qf]:^^F^''d'q  .  (£>.25) 
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Replacing  W  q  with  —  y  and  rearranging  terms  gives 

-  ef  /ln(p  -  qfr^^da,  +  f  j\n{p  -  qfd^U]T^ 

-  \ei7  j  ln(p  -  qfd^^-i^F}"  .  {D.26) 

We  can  now  make  use  of  Eq.(D.13)  and  substitute  the  previous  expression  into 
Eq.(£).9)  to  obtain 

+  f  jlnip-qfd^Ulr^ 
-'^f^Hp-qfd^X.F'^'']-  iD.27) 

At  this  stage,  it  is  possible  to  partially  integrate  the  first  volume  integral  excluding 
the  point  p  =  q  and  therefore 

+    /      P5(p  -  q)-']r^d\  .  (D.2S) 

Jp:^q 

Notice  that  the  last  integral  in  Eq.(I?.28)  vanishes  as  can  be  seen  by  looking  at 
Eq.(D.20)  and  the  fact  that  we  are  excluding  the  point  p  =  g  in  the  integration. 
Using  Eq.(i).16)  once  again,  and  replacing  dq^  by  —d'^  we  obtain, 

/    U]T}d^q      rd/r^da.     ^    r  r^da. 
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Subtituting  this  expression  into  Eq.(D.27)  yields  the  result 


(27r)2  7     {p-qy         (27r)2        f  {p  -  qf 


:^\2  jjn^-nX 


-\fj\n{p-qfd^Ulv' 

+  If  JHp  -  q)^d^2,^^"]  ■  (D.SO) 

Since  S{qi)  and  S{q2)  appear  in  Eq.(D.30)  in  F-^" ,  we  might  expect  singu- 
larities at  (q±  kf  =  n?  =  I]  that  is  at  g^  ±  2^  ■  ^  +  ^^  =  1.  At  this  point  it 
is  convenient  to  make  a  Wick  rotation. ^^  This  rotation  causes  the  real  axis  of  the 
fourth  component  of  the  momentum  to  He  along  the  imaginary  axis.  This  can  be 
done  by  making  the  substitution, 

p^  -^  ip^;  d'^p  =  dp^dp^dp'^dp^  ->  idp^dp^dp^dp'^  .  {D.Zl) 

Making  this  transformation, it  is  possible  to  write  the  condition  at  which  we  can 
expect  a  singularity  as 


r±2q-k  +  k'  =  -sT  2^Jsi-k^)C  -  (-P)  =  1  .  (D.32) 

The  quantity  (  is  the  cosine  of  the  angle  between  the  4- vectors  q  and  k.  where 
~s  —  qfj,q^  =  q^.  If  s  and  {-k"^)  are  always  positive,  the  minimum  value  of  the 
left-hand  side  occurs  at  ±(  =  -1.  Thus,  we  conclude  that 

s  ±  2^s{-P)C  >  s  -  2y^(^+  (-P)  =  [^^  -  ^/IP]2  >  0  , 
and  therefore 


-~sT^yJs(-k^)-(-k')<0<l  .  (D.33) 

Thus,  the  Wick  rotation  and  limitation  to  positive  values  of  s  and  {-k'^)  eliminates 
the  singularities. 
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The  surfaces  over  which  we  integrate  in  Eq.(£).30)  exclude  the  point  (p  =  q) 
by  an  infinitesimal  sphere.  Therefore,  we  can  use  an  infinitesimal  sphere  of  radius 
e'  and  let  this  radius  approach  zero  in  the  end.  Since  da'^  is  the  outward  drawn 
normal,  it  points  toward  the  singularity  at  p  =  ^, 


Ja"  =  _(g''      P'^^'^dQ  .  (D.34; 

Notice  that  all  the  surface  integrals  in  Eq.(I?.30),  except  the  second,  have  the  form 

h^  i  f{\q-p?)9uda^ 


|g-p|  =  e' 


=  -9uf{e)e'^  jiq"  -  p')dQ.  =  0  .  {DM) 


The  second  integral,  however,  is 


r>.j;),,^,,.g./fa-..)r;<iQ 


J\q-p\=e'  (p-q?  """     "  ^  I      {p-q? 

■^Iff— i)l  =  e' 


gr-p|  =  e' 


2r\gi,g2)c^^ 

(p  -  qY 


-2/  1(51,^2^^^.  (r>.36) 

J\q-r\  =  e' 

In  the  limit  5  ^  p,  (  i.e.  e'  — >  0), 

h  =  2r\pi,p2)  /^^  =  2(27r2)r\pi,p2)  .  (i?.37) 

When  these  integrals  are  subtituted  in  Eq.(D.30),  the  vertex  function  V}  on  the 
left-hand  side  is  canceled  leading  to  the  equation, 

0  =  7^  +  ^-1      /^i!£^  +  a./J^^ 

-      (2^)H    L  {p-qf         tiP-qf 

+  l^<f  ln(p  -  qfd^fUlr^  +\ef  i  HP  -  qfF^^'da, 

-  lef  ^Inip-gyd/x^F"^"^  ,  (D.38) 
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where  the  integration  is  over  the  4-dimensional  outer  surface  $R  that  bounds  the 
voliime  in  q-space.  Using  the  identity, 

5''ln(p--,-)^=2^|f|^,  iD.39) 

it  is  possible  to  rewrite  Eq.(Z).38)  in  the  form, 

0  =  7^  +  ^-1    i^/r'^-._2/^V-^^)^"^ 


2    r„  (p-^)2  2     J^  (p-^)^ 


If  we  take  the  surface  3?  to  be  a  large  sphere  of  radius  ^/s  where 

iP  -  q?  =  (p^  -  qt^W  -  qn  -^  q^q"  =  -^ ,  (dai) 

it  follows  that  on  the  sphere  then 

da^  =  s^'^  ^  =s  qudClg  (D.42) 

where  the  symbol  dflg  represents  a  differential  solid  angle  in  q-space.  Substituting 
these  results  in  Eq.(£).40)  gives 

0  =  1^  +  7^1     j  (-srHd/T^)sq,dil  +  2  ^  s-^r^f'sq^dng 

J«£(-*)"'(-/)3/l„£*''ffl,}.  (D.43) 


1 


CanceUing  all  possible  factors  involving  s  and  s    ^  simplifies  Eq.(Z).43)  and  yields 
the  result. 


0  =  7^  +  777^1  -  i  q^{d/T^)dQg  -  2  /  T^dQ^ 
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Notice  that  on  the  integration  surface,  T^  is  constant,  and  therefore  all  terms  in- 
volving derivatives  of  T'^  are  zero.  With  this  simplification  Eq.(D.44)  becomes 

From  Eq.(D.3)  we  see  that  F^"  involves  S  and  T^  where,  by  Eq.(2.37), 

5;-^  =-A  +  B^.  (D.46) 

Since  the  surface  in  the  integral  in  Eq.{DA5)  must  be  pushed  to  infinity,  we  can 
use  the  asymptotic  form  of  5,  namely 

Sr^  -  5/;      5  ~  ^~^B-^  =  }iB-^p-^  .  {DAI) 

From  Eq.(D.3),  it  follows  that 

Z^''(9l,?2)  ~   -  '^iiV\quq2)i2i:'  ,  (DAS) 


and  therefore 


and 

l^F^"  -  -  ^lJiT\q„q2y,Y  .  (D.bO) 

Notice  also  that  in  the  asymptotic  region,  /i  ~  ^2  ~  ^and  therefore  gf  ~  gl  ~  (--s)- 
Substituting  these  forms  into  the  previous  equations  gives, 

q.F^-'  r^(s-')B~'^\q,q)  (D.51) 

and 

l.F^'  -  -  is~')B-'7^yT\q,q)^Y  ■  (D.52) 
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We  assume  that  the  vertex  function  has  the  form, 

L\qi,q2)  =  vx^  +  r\{quq2)  ,  {D.53) 

where    £^(51,52)  — *■  0  as    s  ^-  00.  With  this  assumption,  Eqs.(D.51)  and  {D.52) 
becomes 

and 

Using  the  properties  of  the  Dirac  matrices  (see  appendix  A),  it  is  possible  to  trans- 
form Eq.(jD.55)  into  the  form 

Substituting  Eqs.(Z).54)  and  (D.56)  into  Eq.(Z).45)  gives 
0  =  7'  +  ^{  -  2(27r2)r\  -  ^eB-Y 

-  le[2B-\27r^)x^v  -  ^es'^ B'^  i  q^id^q\  .     (D.57) 
In  obtaining  Eq.(D.58)  we  have  used  the  fact  that 

*  d^q  =  27r2  .  (D.58) 

The  remaining  integral  is  standard  and  is  found  to  be 

/  q^id^q  =  -i^r^  V  •  (I?-59) 

Subtituting  this  result  in  Eq.(D.57)  and  solving  the  equation  for  the  function  v 

gives 


or  since      5  -^  1,  it  follows  that,  to  order  0(e), 


3 

w~l--e.  {DM) 

8 


This  result  is  in  agreement  with  the  findings  of  Yock 
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APPENDIX  E 
HYPERGEOMETRIC  SERIES  AND  ASSOCIATED  LEGENDRE  FUNCTIONS 

E-1    Hypergeometric  series 

If  a  homogeneous  linear  differential  equation  of  the  second  order  has  at  most 
three  singularities,  we  may  assume  that  these  are  at  0,oo,  and  1.  If  all  of  these 
singularities  are  regular,  then  the  equation  can  be  reduced  to  the  form, 

z{l  -  z)—  +  [c  -  {a  +  b  +  l)z]-—  -  abu  =  0  ,  (E.l) 

where  a,  b,  and  c  are  independent  of  z.  This  is  the  hypergeometric  equation.  The 
parameters  a,  6,  and  c  are,  in  general,  complex  numbers. 

Let  us  define 

(a)„  =  r(a  +  n)/r(a)  ,  {E.2) 

i.e., 

(a)o  =  l,  ia)n  =  a{a  +  l)---{a  +  n-l)     n  =  1,2,3,....  {E.3) 

If  c  7^  0,-1, -2,...,  then 

ui  =  5](°)"(^)"  ^"/[(c)n  nl]  =  F{a,  b;  c;  z)  {E.A) 

is  a  solution  of  Eq.(J5'.l)  which  is  regular  at  5r  =  0.  If  c  =  -n  where  n  =  0, 1, 2, . . . , 
then 

oo 
«i  =  z"+  ^(a  +  „  +  1)^(6  +  n  +  1)^  ^-/[(n  +  2)^  m!] 
m=0 
=  z"+^F(a  +  n  +  l,6  +  n  +  l;n  +  2;z)  (£.5) 
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is  a  solution.  The  function  F{a,  b;  c;  z)  is  called  the  hypergeometric  series  of  variable 
z  with  parameters  a,  6,  and  c.  If  c  =  -m,  {m  =  0,1,2,...)  Eq.(£;.4)  becomes 
meaningless.  If  a  =  -n  or  6  =  -n  where  n  =  0, 1,2, ... ,  and  if  c  =  -m  where 
m  =  n,n  +  l,n  +  2, ...,  then  we  define 

F{~n,b;-m;z)  =  Er=o(-'^)r(^)r^7[(-^)r'^!] 

(E.6) 
F{a, -n; -m;  z)  =  Er=o(«)r-(-")r^7[(-"^)r"^!] 

Since  Eqs.{E.5)  and  (E.6)  give  solutions  of  Eq.{E.l),  we  see  that  the  hypergeometric 
equation  has  a  solution  which  is  a  polynomial  in  z  whenever  —a  or  —6  is  a  non- 
negative  integer.(If  a  =  -?ri  ov  b  —  -m  and  c  =  -n,  where  n  =  0,  1,  ...  ,and 
m  =  n  +  1,  n  +  2,  . . . ,  the  series  in  Eq.(£;.5)  terminates.) 

If  a  and  6  are  different  from  0,  -1,  -2,  . . . ,  then  the  hypergeometric  series 
given  in  Eq.(£;.4)  (or  (£'.5),  in  the  case  c  =  -n)  converges  absolutely  for  all  values 
of  \z\  <  1.  For  1^1  =  1  we  have: 

absolute  convergence  for    |^|  =  1    if  Re{a  +  6  —  c)  <  0, 

conditional  convergence  for    |z|  =  1     z  /  1    if  0  <  Re(a  +  b  —  c)  <  1, 

divergence  if    \z\  —  1    and    1  <  Re{a  +  6  —  c). 

The  n*''-derivative  of  equation  (E.4)  is  given  by  the  relation, 

—  F(a,  b;  c;  z)  =  ^     "^  ^"  F{a  +  n,b  +  n;c  + n;  z)  .  {E.7) 

The  hypergeometric  series  satisfies  the  recurrence  relation, 

(c  -  a)(c  -  b)zF{a,  b;c  +  l-z)  =  c  [(2c  -  a  -  b  -  l)z  -  c  +  1]  F{a,  b;  c;  z) 

+  c{c  -  1){1  -  z)F(a,b;  c  -  1;  z)  .        (E.S) 
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E-2   Associated  Legendre  Functions 

The  Legendre  functions  are  solutions  of  the  differential  equation, 

where  //,  i/,  and  z  are  unrestricted. 

1 
Under  the  substitution,  w  =  [1  —  z'^)'i^v,  Eq.(£'.9)  becomes 

(1  -  z^)—  -  2(^  +  l)z—  +  (z/  -  n){v  +  ^  +  l)y  =  0  .  {E.IQ) 

With  (  =  I  —  |z  as  the  independent  variable,  this  differential  equation  becomes 
d^X)  dv 

This  is  the  same  as  Eq.{E.l)  if  we  set  a  =  /x  —  i/,   6  =  /u  +  i^  +  1,  and  c  =  pL  +  1. 
Hence  it  follows  that  the  function, 


V 


^"W  =  f(T^(f^)'  n-.-  +  i;i-.;|-J^), 


{E.12) 


is  a  solution  of  Eq.(E.lO)  provided  |1  -  2;|  <  2.  The  solutions  P^(^)  are  called 
the  Associated  Legendre  Functions  of  the  first  kind.  If  we  set  (  —  z'^,  Eq.(E.lO) 
becomes 

4C(l-C)^  +  [2-(4C  +  6)C]^-(/.-t^)(/.  +  z^  +  l>  =  0  (£;.13) 

which  is  also  of  the  hypergeometric  type  with  a  —  |(/i  +  u  +  1),  b  —  |(/i  —  i/),  and 
c  =  |.  Hence  Eq.(^.l)  has  a  solution  of  the  form, 

w  =  Q':;{z)  =  e'f^^2-''-'V^^^^^-^^^^z-'->'-\z^  -  1)^'^ 

X  F(|r/  +  Im  + 1,  Ia'  +  1'^  +  ^;  z^  +  |;  z~^)  ■       (E.u) 
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The  solutions  Q'^{z)  are  called  the  Associated  Legendre  Functions  of  the  second 
kind.  Legendre's  differential  equation  remains  unchanged  if  fi  is  replaced  by  — ^,  z 
hj  —z,  and  u  hy  —u  —  1.  Therefore 


P±'^(±z),     Qfi^z),     P-^i(±^),     and     QS-i(±^) 


(^.15) 


are  solutions  of  Eq.(£'.l).  Using  various  transformation  formulas  of  the  hyperge- 
ometric  series^^  it  is  possible  to  write  these  functions  in  several  ways.  These  are 
presented  in  Tables  E.l  and  E.2.  The  Wronskian  of  the  functions  P^{z)  and  Q^{z) 
is  given  by  the  expression, 


The  Legendre  functions  satistify  the  following  relations: 


{E.16) 


P!:{z)^p>t^_,{z)- 

^'^'^  =  r(l|!!-ii)^'-"'^"^  +  (2/7r)e-^'^'^  sin(^^)Q^(.)  ; 
Pu'{^)  =  rl j  "!!!!!!  {Pi:{^)  -  (2^)6—"  sin(/.vr)g^(.)}    ; 


p;^'{^)  = 


sin[7r(^-^)]{Q(:(^)-QVi(^)} 


TV  cos(i/7r)r(l  +  ^  +  /i) 
(5^(^)sin[7r(i/  +  /i)]  -  g::^_i(z)sin[7r(i/  -  ^)]  =  7re'>'^  cos(i/7r)P,''(^)  ; 


Q(;(z)sin(7r^)  =  i7re'>'^ 


gl,_i(-~)  -  Qtiz)  =  e*-''  cos(t/7r)r(l  +  //  +  i^)r(//  -  tv)p-''(z) 


{E.17) 

(E.IS) 

{E.19) 

(£.20) 
(E.21) 
(£;.22) 
(£.23) 
(£.24) 
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Table  E-1    Expansions  for  P/'(z) 


P^i^)  =  r(^  (^)  ''  (.  +  lYFi-u,  -.-,;!-.;  i^)  (^.26) 

p'(')  =  Yxi~;r)^''~'^^~''''^^^  (^-27) 


p.-'  ^ 


'.''(^)  =  f(rr7;)(^'  - 1)"^'[^  +  V^^^r-' 


1 


X  ^(-^^  -  Ia^,  I  -  //,  1  -  2/.;  ^5^^)  (£;.29) 

X  F{-u  -  ^,  I  -  /.,  1  -  2^;  -;_^(gg^)  (E.30) 


P.'i^)  =  ~y-^ {z>^—\z^  -  1)-^" 


3.  ^-2> 


X  P(-|/i  _  ij.,  i  -  li.  -  1^,  l-u-  z-'')  .  (E.31) 
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Table  E-2    Expansions  for  e^2'^''Q{;(z) 


xF{l  +  fx  +  u,l  +  u;2u  +  2;^J  (E.32) 

X  ^(1  +  1^  -  |a^,  I  +  1^  +  |)«,  ^  +  I;  T~^)        {E.33) 


—  X  —  77  1/ 


X  F(  J  +  |i/  +  i/i,  -ii/  +  1^,  1  +  ^;  1  -  ^2)  (_£;_36) 
Several  values  of  /i  and  v  are  of  interest  and  are  presented  below. 

Q3(^)  =  z(2/7r)2(22  _  ^y\^^  ^  y^2_l]-.'-i  (^_3g) 

p-'(z)  =  (2/7r)^^^^^i^|[^  +  y^^l]-+i  -  [^  +  yi^^l-^-ij  (£;.39) 


qH^)-- 

z2-l)-i 
2u  +  l 
2-u 

QH^)-- 

O.  .      1      1 
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V^^-l]"""^  (^.40) 

-(^'-1)^"  (E.41) 

The  contiguous  Legendre  functions  satisfy  the  following  relations: 

Pi:^\z)  +  2(/i  +  i)z{z^  -  i)-i/2  pi:+\z)  =  {v-  ^)iu  -1,  +  i)Pi:{z)  ■  {eaz) 

=  {2u  +  l)^Jz^-lPI^+\z)  ■  (EA5) 

P'-ii^)  -  Pi:+i{z)  -  -(2^  +  lWz^-lPr\^)  ;  (^.46) 


P'-iU)  -  'P'(^)  --{^-f^  +  l)^/^^^Pr\z)  ■  {EA7) 


zPI^{z)  -  Pi:^,{z)  =  -(u  +  i,)Vz^-lPr\z)  ;  {EA8) 

{v  -  f^)zPi;(z)  -{u  +  ^i)P:_,{z)  =  ^/z^-lPi:+\z)  ■  (E.49) 

{u-^x  +  1)P,VW  -  (^^  -  /^  +  1>P.^(^)  =  \/ z^  -  I  Pi:+\z)  ■  (E.50) 

(^2  _  i)^Mf)  =(._;,  +  l)P;+,(.)  -  (a.  +  i^>P.-(.~)  ; 

^vzPi:{z)-{^i  +  v)Pi:_,{z).  (^.51) 

Using  the  expressions  given  in  Tables  E-1  and  E-2  it  is  possible  to  find  the  behavior 
of  the  Legendre  functions  of  the  first  and  second  kind  near  the  singularities,  z  =  1 
and  2;  =  oo.  The  results  are  presented  in  Table  E-3. 
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Table  E-3    Behavior  of  Pi^{z)  and  Qt{z)  at  the  Singularities. 

Behavior  at  2  =  1 

P!^{^)  -  f(r^^^'  ~  ^^"'^         A'  7^  1,  2,  3, ... ,  {E.b2) 

^^^'^  ^  V{u-m^l)^'  -  ^^"'"^         m  =  0,  1,  2, . . . ,  (£;.53) 

e"  =  "''g:;(^)~2-i+  =  ''r(^)(z-l)-^'^  i?e/x  >  0  (E.54) 

Behavior  at  2;  =  00 


(£.55) 
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